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Abstract 

Unmanned  aerial  vehieles  (UAVs)  have  onboard  sensors  that  eontinuously  reeord 
weather  data  during  their  missions.  This  information  is  extremely  valuable  to  both  the 
meteorologieal  and  UAV  eommunities  with  numerous  potential  benefits,  whieh  inelude 
improved  weather  foreeast  produets  and  additional  weather  intelligenee  for  military 
planners.  The  value  of  any  dataset  is  direetly  related  to  its  aeeuraey  and  this  researeh 
determined  the  aeeuraey  of  weather  data  obtained  from  a  partieular  UAV,  the  RQ-4A 
Global  Hawk.  This  was  aeeomplished  through  statistieal  analysis  and  eomparisons  with 
upper-air  data  and  Atmospherie  Slant  Path  Analysis  Model  (ASPAM)  profiles  of  the 
atmosphere.  Reeommendations  are  provided  for  the  use  of  this  valuable  environmental 
intelligenee  souree  to  multiple  user  eommunities. 

Similar  sensors  exist  on  eommereial  aireraft  using  the  Aireraft  Communieations 
Addressing  and  Reporting  System  (AGARS).  Data  from  ACARS-equipped  aireraft  are 
eompiled  and  quality  eontrolled  by  the  National  Weather  Serviee  Foreeast  Systems 
Laboratory,  then  proeessed  and  made  available  to  numerous  ageneies.  Personnel  use  the 
information  not  only  to  enhanee  their  foreeast  produets  but  also  as  a  data  souree  for 
ingest  into  numerieal  weather  predietion  models.  AGARS  data  are  more  spatially  and 
temporally  available  than  rawinsonde  data,  thus  potentially  having  a  more  signifieant 
impaet  on  upper-air  analysis  models.  Foreeasters  also  use  this  near-real-time  data  to 
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enhance  their  products  such  as  weather  warnings  and  advisories.  ACARS  information  is 
a  proven  asset  to  the  weather  community  as  well  as  mission  planners. 

Methods  analogous  to  the  implementation  and  quality  control  of  ACARS  data  can 
also  be  applied  to  the  information  obtained  from  UAVs  since  its  accuracy  is 
demonstrated  in  this  research.  This  study  illustrates  the  utility  of  the  UAVs  to  create  on- 
demand  upper-air  soundings  for  any  location  worldwide,  whereas  ACARS  data  are 
limited  to  commercial  aircraft  routes.  The  use  of  this  data  can  greatly  enhance  forecast 
models  and  products  especially  in  data-sparse  regions,  as  well  as  provide  a  weather 
reconnaissance  capability  to  military  planners. 
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VERIFICATION  OF  METEOROFOGICAF  DATA  REPORTS 

FROM  THE  RQ-4A  GFOBAE  HAWK  UNMANNED  AERIAL  VEHICLE 

1.  Introduction 

Unmanned  aerial  vehieles  (UAVs)  ean  provide  eontinuous  weather  data  during  all 
aspeets  of  their  flight  mission.  This  information,  if  made  available,  would  be  an  asset  not 
only  to  weather  foreeasters  but  also  UAV  operators  and  mission  planners.  An  inerease  in 
the  amount  of  environmental  intelligenee  would  improve  weather  products  and  aid  UAV 
mission  planners  and  operators  in  their  decision-making  processes.  The  value  of  any 
dataset  is  directly  proportional  to  its  accuracy.  Irrelevant  or  inaccurate  information 
would  only  amplify  errors  already  inherent  in  any  process.  Therefore,  efforts  must  be 
made  to  ensure  the  validity  of  available  data  from  UAVs. 

Similar  sensors  exist  on  some  commercial  aircraft,  which  transmit  weather 
information  along  with  other  flight  data  by  a  method  called  Aircraft  Communication 
Addressing  and  Reporting  System  (ACARS).  The  environmental  intelligence  obtained 
from  these  aircraft  has  been  used  successfully  for  the  past  decade  to  improve  numerical 
weather  prediction  (NWP)  models,  enhance  forecast  products,  and  assist  in  flight 
planning  (Moninger  et  al.  2003).  Pressure,  temperature,  wind  speed,  and  wind  direction 
are  reported  by  all  ACARS  equipped  aircraft  and  some  more  advanced  aircraft  provide 
moisture  and  turbulence  data  as  well.  There  are  numerous  other  uses  for  these  data  and 
studies  continue  to  show  their  value.  Before  this  information  was  made  available  to  the 
weather  community  and  other  agencies,  it  was  examined  to  determine  its  accuracy.  It 
was  shown  to  be  a  valid  source  of  environmental  information  and  then  provided  for  use. 


This  thesis  follows  the  style  of  the  American  Meteorological  Society  (AMS  2002). 
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Quality  control  is  a  key  ingredient  in  the  proeessing  of  the  data  and  ensures  the  eontinued 
aeeuracy  of  the  information.  Thus,  the  value  of  produets  ereated  with  this  information  is 
maintained. 

UAV  meteorological  data  may  be  used  and  validated  in  mueh  the  same  way  as 
ACARS  weather  data.  This  study  briefly  deseribes  the  ACARS  proeess  and  presents 
details  on  weather  observations  provided  by  the  system.  Some  of  the  uses  and  benefits  of 
the  data  are  also  explained.  Then  a  short  diseussion  on  eurrent  quality  eontrol  methods 
will  illustrate  the  validity  of  these  data.  A  eomparison  is  then  made  between  ACARS  and 
UAV  data,  showing  how  environmental  intelligenee  available  on  UAVs  ean  be  exploited 
to  improve  weather  forecast  products  and  support.  The  methodology  of  this  researeh 
demonstrates  the  verifieation  of  UAV  weather  data  and  reeommendations  state  that  this 
information  ean  and  should  be  used  in  ways  comparable  to  the  methods  applied  to 
eommereial  aireraft  ACARS  data. 

a.  Problem  statement 

Aceurate  environmental  information  from  UAV  platforms  is  important  to  the 
meteorologieal  eommunity  as  well  as  mission  planners  and  operators.  Military  UAVs 
have  the  ability  to  provide  weather  intelligence  in  data-sparse  regions  of  eombat 
operations,  while  ACARS  equipped  eivilian  aireraft  are  restrieted  from  these  areas.  The 
RQ-4A  Global  Hawk  (Figure  1)  eurrently  reports  weather  data  but  this  information  has 
never  been  validated.  The  following  statement  is  from  the  88**'  Weather  Squadron’s 
(2002)  thesis  topie  proposal,  whieh  summarizes  the  seope  of  this  problem: 
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“In  an  August  27,  2002,  article  in  Federal  Computer  Week,  the  CSAF,  General 
John  P.  Jumper,  said  the  Air  Force  should  begin  to  think  of  the  YF-22  fighter  jet  as 
also  being  an  intelligence,  surveillance  and  reconnaissance  aircraft  and  ‘get  away 
from  being  platform-centric  ’  (Caterinicchia  2002).  In  this  same  spirit,  specialized 
surveillance  platforms  such  as  the  RQ-4A  Global  Hawk  could  theoretically  be 
exploited  to  obtain  environmental/meteorological  information  over  data-sparse 
areas.  These  data  could  then  be  used  by  DoD  to  initialize  meteorological  models. 
There  are  a  number  of  questions  however  that  must  be  answered  before  AFWA  can 
begin  incorporating  meteorological  data  from  platforms  such  as  RQ-4A  Global 
Hawk  into  its  models.  Among  these  questions  is  how  accurate  are  the  data?  ” 


Figure  1.  The  RQ-4A  Global  Hawk  (Northrop  Grumman  2001). 


This  researeh  validates  Global  Hawk  weather  data  using  statistieal  comparisons 
with  other  available  environmental  information.  This  thesis  provides  recommendations 
for  the  use  of  this  information  once  its  accuracy  has  been  shown.  The  Atmospheric  Slant 
Path  Analysis  Model  (ASP AM)  is  the  primary  tool  used  to  verify  the  Global  Hawk  data. 
ASP  AM  is  a  model  used  by  the  Air  Force  Combat  Climatology  Center  (AFCCC)  and 
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ingests  all  available  weather  information  to  ereate  a  “best  guess”  vertieal  profile  of  the 
atmosphere  by  optimal  interpolation  teehniques.  The  data  used  by  ASPAM  inelude 
surfaee  and  upper-air  observations,  ACARS  and  other  aireraft  reports,  satellite  soundings, 
and  information  from  other  meteorologieal  databases,  primarily  the  Navy  Operational 
Global  Atmospherie  Predietion  System  (NOGAPS)  model  analyses. 

b.  Objectives 

The  primary  objeetive  of  this  researeh  is  to  validate  Global  Hawk  weather  data 
and  justify  the  benefits  of  using  sueh  information  in  weather  predietion  models,  weather 
foreeasting,  and  as  weather  intelligenee  for  mission  planning.  Speeifie  steps  needed  to 
aehieve  this  overall  goal  are  to; 

1)  obtain  and  format  environmental  data  from  several  Global  Hawk  flights; 

2)  obtain  and  format  vertieal  atmospherie  profiles  from  ASPAM; 

3)  obtain  and  format  rawinsonde  observations  (RAOBs)  taken  from  Edwards  Air 
Foree  Base  (AFB)  near  the  time  of  the  Global  Hawk  flights; 

4)  eompare  Global  Hawk  weather  data  with  ASPAM  and  RAOB  observations; 

5)  determine  the  validity  of  Global  Hawk  data  based  on  these  eomparisons  and 
statistieal  analyses; 

6)  and,  make  eonelusions  and  reeommendations  to  users  based  on  the  validity  of 
the  Global  Hawk  data. 
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c.  Approach 

1)  Data  ACQUISITION 

Weather  data  were  obtained  from  several  Global  Hawk  flights  in  order  to  perform 
statistieal  analyses.  The  Global  Hawk  System  Program  Offiee  (SPO)  at  Wright-Patterson 
Air  Foree  Base  (WPAFB)  funded  travel  to  Edwards  AFB,  California  to  gather  data  from 
Global  Hawk  flight  tests.  Unfortunately,  flight  tests  did  not  oeeur  during  this  temporary 
duty  (TDY)  but  data  were  aequired  from  missions  performed  earlier  in  the  year. 
Information  was  also  available  from  operational  missions  but  the  use  of  that  data  was 
impraetieal  beeause  of  the  sensitive  nature  of  the  operations.  Data  obtained  from 
elassified  operations  and  missions  flown  in  sensitive  areas  are  elassified  and  diffieult  to 
obtain  and  use.  Weather  data  were  also  obtained  from  ASPAM  analyses  and  Edwards 
AEB  RAOBs  to  eompare  with  the  Global  Hawk  observations.  The  times  and  loeations  of 
the  ASPAM  profiles  were  provided  at  the  medians  of  the  Global  Hawk  aseent  and 
deseent  times  and  loeations,  so  that  the  data  were  as  representative  as  possible.  Edwards 
RAOB  data  were  provided  at  the  elosest  possible  dates  and  times  to  the  Global  Hawk 
missions  but  were  unavailable  at  the  speeifie  times  of  certain  missions. 

2)  Data  FORMATTING 

The  Global  Hawk  weather  information  was  separated  and  decoded  from  the  data 
stream.  This  was  initially  accomplished  during  the  Australian  deployment  of  the  Global 
Hawk  in  2001  and  similar  methods  were  investigated  during  this  research.  Data 
processors  of  the  452"‘*  Elight  Test  Squadron  (EETS)  decoded  and  compiled  data  from 
five  separate  missions  and  provided  the  information  for  this  study  on  CD-ROM  media. 
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The  Global  Hawk,  ASPAM,  and  RAOB  data  were  then  put  into  the  same  format  and  unit 
eonversions  were  performed. 

3)  Data  VERIFICATION 

All  of  the  formatted  data  were  entered  into  statistical  software  packages  and 
several  tests  were  performed.  The  intent  of  these  tests  was  to  show  the  accuracy  of  the 
Global  Hawk  weather  data  when  compared  with  actual  environmental  observations.  It 
was  assumed  that  the  ASPAM  and  RAOB  data  represented  actual  ground-truth 
conditions,  although  ASPAM  interpolates  available  observations  to  a  particular  location 
and  time,  and  RAOB  data  were  not  available  at  the  exact  time  and  location  as  the  Global 
Hawk  information.  The  statistical  measures  used  include  simple  linear  regression, 
analysis  of  variance  (ANOVA),  coefficient  of  determination  (R  ),  mean  square  error 
(MSB),  root  mean  square  error  (RMSE),  and  hypothesis  testing  using  the  t-test. 
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2,  Literature  Review 


a.  ACARS  background 

Weather  observations  are  reported  by  eommereial  aireraft  in  the  United  States 
(US)  through  the  Aircraft  Communications  Addressing  and  Reporting  System  (ACARS) 
which  is  managed  by  Aeronautical  Radio,  Inc.  (ARINC).  The  data  are  compiled, 
processed,  quality  controlled,  and  made  available  by  the  National  Oceanic  and 
Atmospheric  Administration’s  (NOAA)  Forecast  Systems  Laboratory  (FSL).  FSL 
collects  over  170,000  meteorological  observations  daily  including  data  from  aircraft 
outside  the  US,  which  are  known  as  aircraft  meteorological  data  relay  (AMDAR)  reports, 
and  displays  them  on  an  interactive  website  (FSL  2003,  http://www.acweb.fsl.noaa.gov/). 
Pressure,  temperature,  wind  speed,  and  wind  direction  are  the  primary  environmental 
variables  reported,  although  some  aircraft  provide  moisture  and  turbulence  data.  The 
spatial  coverage  of  the  reports  is  shown  in  Figure  2  for  the  world  and  Figure  3  for  the 
continental  United  States  (CONUS)  (Moninger  et  al.  2003). 

The  data  are  currently  accessible  only  by  National  Weather  Service  Forecast 
Offices  (NWSFO)  and  other  government  agencies.  Private  organizations  and  research 
groups  may  obtain  the  information  and  are  encouraged  to  use  it  in  their  research  and 
operations  as  long  as  certain  criteria  are  met  (FSL  2003, 
http://acweb.fsl.noaa.gov/FAO.html). 
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All  dato  for  27-Mar-02  0000  to  27-Mar-02  2359  UTC 

130070  reports  shown,  between  —500.  and  43013.  ft 


Figure  2.  A  24-h  worldwide  ACARS/AMDAR  plot  for  27  Mar  2002.  ACARS  reports 
are  in  red  and  AMDAR  reports  are  in  blue  (Moninger  et  al.  2003). 
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Figure  3.  ACARS  coverage  over  the  CONUS  for  27  Mar  2002.  Reports  are  between 
25,000  and  45,000  ft  and  are  color  coded  by  altitude  (Moninger  et  al.  2003). 
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b.  ACARS  data  use 


ACARS  and  AMDAR  information  have  numerous  uses  and  benefits  to  the 
meteorological  community  such  as  improved  NWP  potential,  enhanced  weather  forecast 
products,  and  an  overall  increase  in  available  upper-air  observational  data.  ACARS 
reports  are  also  more  spatially  and  temporally  available  than  current  rawinsonde 
information.  Moninger  et  al.  (2003)  state  the  “data  are  the  only  in  situ  source  of  upper-air 
data  at  non-synoptic  times,  and  therefore  are  highly  valuable.”  Some  of  the  benefits  of 
this  data  are  illustrated  in  the  following  case  studies: 

1)  Thunderstorms 

Mamrosh  (1998)  showed  the  value  of  ACARS  data  in  forecasting  a  convective 
event  by  the  Chicago  NWS  office.  ACARS  environmental  information  was  used  to 
determine  the  amount  of  atmospheric  instability  and  wind  shear  during  the  day  and 
accurately  predict  the  timing  and  type  of  severe  weather  to  occur. 

Figure  4  is  a  sounding  created  from  ACARS  data  at  1515  Universal  Time  Code 
(UTC),  which  showed  an  inversion  indicating  a  stable  layer  in  the  lower  atmosphere. 

The  convective  temperature  was  established  with  the  data  and  the  forecasters  believed  it 
was  possible  to  reach  this  temperature  due  to  strong  warm  air  advection  in  the  region. 
Also  evident  on  the  sounding  was  a  low-level  jet  of  60  knots,  which  could  mix  down  to 
the  surface  once  the  inversion  broke.  Figure  5  is  an  ACARS  sounding  about  one  hour 
later  which  shows  the  inversion  had  lifted  and  weakened  and  the  low-level  winds  were 
still  significant  at  40  to  60  knots. 
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Descent  sounding  from  234°  into  Chicago/Ohare,  iL  (ORD) 


Figure  4.  ACARS  sounding  from  O’Hare  at  1515  UTC  (Mamrosh  1998). 


Ascent  sounding  toward  82°  from  Chicago/Ohare,  iL  (ORD) 


Figure  5.  ACARS  sounding  from  O’Hare  at  1623  UTC  (Mamrosh  1998). 
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A  sufficient  amount  of  instability  to  initiate  thunderstorm  aetivity  was  evident  on 
a  sounding  about  an  hour  later  (Figure  6)  and  a  Severe  Thunderstorm  Watch  was  issued 
by  the  Storm  Prediction  Center  (SPC).  Thunderstorms  developed  in  the  area  shortly 
thereafter  and  several  severe  weather  events  followed  such  as  winds  of  50  to  60  knots  and 
golfball-sized  (1.75  inch)  hail  (Mamrosh  1998).  Mamrosh  (1998)  stated,  “the  ACARS 
data  was  a  very  useful  supplement  to  the  usual  tools  found  in  an  NWS  office,  for  there  is 
no  other  system  in  place  to  get  very  frequent  soundings  of  both  wind  and  temperature.” 


Descent  sounding  from  281°  into  Chicago/Ohare,  iL  (ORD) 


Figure  6.  ACARS  sounding  from  O’Hare  at  1738  UTC  (Mamrosh  1998). 
2)  Turbulence 


Reports  from  aircraft  using  ACARS  were  extremely  benefieial  in  verifying 
turbulence  for  the  Center  Weather  Service  Unit  (CWSU)  at  Seattle  (Mamrosh  et  al. 
2001).  The  polar  jet  stream  was  loeated  across  the  Pacific  Northwest  with  maximum 


11 


winds  in  excess  of  175  knots  through  Oregon  (Figure  7).  Rawinsonde  data  at  the  level  of 
the  jet  stream  were  unavailable  because  the  high  winds  transport  the  instrument  far  from 
the  ground  station  and  the  signal  is  lost.  Therefore,  ACARS  data  were  the  only  existing 
information  and  were  used  to  verify  the  placement  of  the  jet  stream  and  determine  areas 
of  moderate  to  severe  turbulence.  Additional  ACARS  reports  also  allowed  the  Seattle 
CWSU  to  track  the  slow  northward  progression  of  the  jet  stream  during  the  day. 


Figure  7.  ACARS  wind  data  between  25,692  and  37,008  ft,  for  the  period  1200  UTC  to 
1459  UTC  January  17,  1998  (Mamrosh  et  al.  2001). 


Some  aircraft  also  report  turbulence  through  varying  measures  and  algorithms 
including  Derived  Equivalent  Vertical  Gust  (DEVG),  Eddy  Dissipation  Rate  (EDR),  and 
measuring  the  vertical  acceleration  of  aircraft  (Moninger  et  al.  2003).  Some  of  this 
information  has  inaccurate  values  and  biases  in  the  algorithms  but  the  data  are  useful 
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nonetheless.  EDR  is  the  recommended  and  most  valuable  measure  of  turbulence  because 


it  relates  to  the  atmosphere  rather  than  the  aircraft.  The  other  two  procedures  rely  on  the 
weight  and  velocity  of  the  aircraft  to  determine  a  value  for  turbulence. 

3)  Winter  WEATHER 

Determining  whether  precipitation  will  be  liquid  or  solid  in  winter  storms  is 
always  a  challenge,  especially  with  minimal  data.  One  such  occasion  where  ACARS 
reports  were  used  to  improve  winter  forecasts  was  during  a  moderate  snow  event  in 
Chicago  on  Christmas  Eve  (Mamrosh  et  al.  2001).  Most  of  the  NWP  models  predicted 
rain  for  the  24**'  and  25**^  of  December,  while  others  hinted  at  only  a  trace  of  snowfall 
during  the  period.  Two  soundings  from  ACARS  data  (Eigures  8  and  9)  were  examined 
and  a  layer  with  temperatures  below  freezing  was  evident  in  the  low  levels.  The 
soundings  also  showed  a  mid-level  layer  with  above  freezing  temperatures  and  that  its 
thickness  was  decreasing.  Wind  data  from  the  soundings  showed  low-level  easterly  flow, 
which  would  advect  dry  air  into  the  region  causing  the  air  to  cool  even  more  as  it 
evaporated.  This  information  indicated  the  possibility  of  freezing  rain  and  snow  in  the 
area  and  forecasts  were  updated  as  such. 

The  precipitation  began  that  evening  as  freezing  rain  and  changed  over  to  snow 
soon  thereafter.  The  snow  continued  overnight  and  ended  as  drizzle  the  following 
evening.  The  amount  of  snowfall  recorded  was  about  2.5  inches. 


13 


Figure  8.  ACARS  sounding  showing  above-freezing  layer  (Mamrosh  et  al.  2001). 


Figure  9.  ACARS  sounding  showing  above-freezing  layer  diminishing  (Mamrosh  et  al. 
2001). 
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Mamrosh  et  al.  (2001)  stated,  “This  was  an  especially  useful  forecast  as  even 
modest  snowfalls  as  this  has  significant  impact  on  a  busy  airport  such  as  O’ Hare  - 
especially  the  day  before  Christmas!”  Without  the  ACARS  reports,  this  forecast  might 
have  been  missed. 

4)  Airline  MISSION  CONTROL 

Another  benefit  of  ACARS  data  is  that  they  are  near  real  time  and  can  assist  flight 
planners  and  mission  controllers  with  weather  related  decisions.  An  example  of  this 
occurred  when  an  Italian  airliner  encountered  unexpected  strong  headwinds  over  the 
Atlantic  Ocean  and  the  crew  unsure  if  the  fuel  supply  would  be  sufficient  to  complete  the 
flight.  The  Miami  CWSU  used  ACARS  wind  data  to  determine  an  altitude  with  lighter 
headwinds  (Figure  10),  the  airliner  was  directed  there  by  controllers,  and  the  flight  was 
completed  without  incident  (Mamrosh  et  al.  2001). 


Figure  10.  ACARS  wind  data  between  22,000  and  39,000  ft  (Mamrosh  et  al.  2001). 
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The  case  studies  above  demonstrate  the  value  of  ACARS  data  in  the  forecasting 
process,  which  is  just  a  hint  of  the  possible  uses  and  benefits  of  the  information.  The 
reports  are  also  used  with  NWP  models  both  as  input  for  enhanced  model  initialization 
and  as  verification  for  model  forecasts. 

c.  ACARS  accuracy  and  quality  control 

ACARS  weather  data  are  only  valuable  if  they  are  accurate  and  checked  for 
errors.  Several  studies  have  shown  the  accuracy  of  ACARS  reports  by  comparing  them 
with  other  sources,  such  as  rawinsondes  and  models  (Moninger  et  al.  2003).  ACARS 
data  were  compared  with  rawinsonde  readings  by  Schwartz  and  Benjamin  (1995)  and  the 
values  were  shown  to  agree.  Further  statistical  analysis  implied  ACARS  reports  were 
actually  more  precise  than  rawinsonde  data.  Benjamin  et  al.  (1999)  determined  wind  and 
temperature  errors  in  the  400-300  millibar  layer  to  be  1.1  m  s'^  and  0.5  °C,  respectively. 
The  low  variability  of  mesoscale  values  was  the  reason  this  atmospheric  layer  was  used. 

Once  the  data  were  shown  to  be  accurate,  they  were  made  available  for  use  in 
forecasting  and  mission  control.  Continuous  quality  checks  must  be  performed  on  the 
data  to  ensure  sustained  accuracy  and  value.  Quality  control  statistics  from  a  system  used 
by  the  U.S.  Navy  is  seen  in  Table  1  (Pauley  2003,  in  preparation).  The  FSL  utilizes  a 
quality  control  system  that  checks  reports  and  notifies  the  airlines  of  possible  equipment 
problems.  When  reports  from  an  aircraft  are  consistently  in  error,  the  airline  is  advised 
so  corrections  can  be  made.  This  relationship  between  the  data  users  and  providers  is 
important  to  guarantee  continued  value  of  ACARS  information. 
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Table  1.  Percentage  of  rejected 
(Moninger  et  al.  2003). 

data  from  U.S. 

Navy  quality  control  methods 

ACARS 

AMDAR 

Duplicates 

2.1% 

0.8% 

Bad  reports 

2.5% 

1.1% 

Bad  temperature 

0.5% 

0.0% 

Bad  winds 

0.8% 

0.4% 

d.  UA  V  background 

Unmanned  aerial  vehieles  report  environmental  information  much  the  same  way 
as  commercial  aircraft  but  the  data  are  not  as  readily  accessible.  Comparisons  can  be 
made  between  UAV  and  ACARS  reports  showing  the  potential  value  of  the  weather 
information.  Unmanned  systems  perform  a  variety  of  missions  at  varying  locations  and 
altitudes  and  for  different  lengths  of  time.  The  environmental  intelligence  would  benefit 
both  the  weather  forecasters  and  the  UAV  controllers  by  providing  much  needed  weather 
information  in  data-sparse  regions,  which  is  where  most  of  the  flights  occur. 


e.  Possible  uses  of  UA  V  data 

Forecasters  could  use  UAV  weather  information  to  improve  their  products  and 
enhance  the  NWP  models.  The  increased  amount  of  information  available  to  weather 
personnel  can  be  utilized  in  ways  comparable  to  ACARS  reports,  as  seen  in  Section  2.b. 
Currently,  the  U.S.  Navy  uses  ACARS  data  in  their  synoptic  scale  model  NOGAPS  and 
their  mesoscale  model  called  the  Coupled  Ocean/ Atmospheric  Mesoscale  Prediction 
System  (COAMPS)  and  are  interested  in  obtaining  UAV  data  as  well  (E.C.  Mozley  2002, 
personal  communication).  The  U.S.  Air  Force  has  also  shown  interest  in  using  UAV  data 
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in  their  mesoscale  model  ealled  the  Weather  Researeh  and  Foreeasting  model  (WRF). 
The  data  eould  also  alert  UAV  controllers  of  possible  turbulence  and  low  temperatures, 
which  impact  the  performance  and  stabilization  of  the  UAV  systems.  Another  advantage 
to  the  operators  is  using  the  data  to  redirect  missions  to  areas  with  weather  more 
favorable  to  the  aircraft  and  its  sensors  and  weapons. 

/  UA  V  accuracy  and  quality  control 

Weather  data  from  UAV  sensors  should  be  checked  for  accuracy  similar  to 
methods  used  with  ACARS  reports.  The  information  should  be  made  available  for  use 
once  the  validity  of  the  data  is  shown.  Then  quality-control  procedures  must  be 
established  and  followed  to  ensure  the  continued  value  of  the  data.  Using  procedures 
which  are  analogous  to  ACARS  methods  is  a  feasible  approach.  The  Fleet  Numerical 
Meteorology  and  Oceanography  Center  (FNMOC)  has  advanced  algorithms  in  place  to 
test  data  used  in  NOGAPS  for  accuracy  and  relevancy.  These  and  other  similar  methods 
could  be  used  with  UAV  weather  data. 
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3,  Data  Collection  and  Review 


a.  Data  collection 

Three  sets  of  data  were  analyzed  and  eompared  in  this  research.  The  primary 
dataset  was  from  the  Global  Hawk,  which  included  the  environmental  parameters  that 
were  evaluated.  Vertical  atmospheric  profiles  from  the  Atmospheric  Slant  Path  Analysis 
Model  (ASPAM)  and  upper-air  rawinsonde  soundings  from  Edwards  Air  Force  Base 
(AFB)  were  used  to  simulate  “ground  truth”  datasets  that  were  compared  with  the  Global 
Hawk  data  observations. 

1)  Global  Hawk 

All  Global  Hawk  datasets  used  in  this  research  were  obtained  within  the  Edwards 
AFB  operating  range  area.  The  positions  varied  from  34.82  to  36.56  degrees  north 
latitude  and  116. 6  to  118. 6  degrees  west  longitude  (Figures  1 1  and  12).  Observations 
were  obtained  from  five  separate  test  flights  occurring  in  2003  on  29  May,  3  June,  26 
June,  8  August,  and  15  August.  All  data  were  from  the  same  aircraft.  Air  Vehicle  7 
(AV7),  which  is  the  newest  airframe  in  the  inventory.  The  intent  was  to  obtain  data  near- 
real-time  while  on  TDY  to  Edwards  AFB  during  late  September  2003.  Flight  testing, 
however,  did  not  occur  during  this  period  due  to  weather  or  to  be  more  precise,  the  lack 
of  weather.  The  452”‘*  FFTS  had  planned  to  perform  crosswind  tests,  which  require  wind 
speeds  of  10  to  15  knots  for  takeoff  and  landing,  but  an  area  of  high  pressure  dominated 
the  region  and  winds  were  unusually  calm  during  the  week.  Therefore,  the  only  available 
data  were  from  previous  missions  and  the  time  at  Edwards  was  spent  working  with  the 
data-processing  section  to  get  the  information  that  was  used  in  this  study.  The  desired 
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parameters  were  determined  and  a  request  was  made  to  the  412*  Test  Wing  Range 
Division  at  Edwards  AFB,  which  provided  the  information  in  comma-separated  text 
format  on  CD-ROM  media. 

Data  points  for  each  mission  were  provided  once  per  second,  which  was  more 
than  adequate  for  this  study.  The  duration  of  the  missions  ranged  from  one  hour  and  45 
minutes  to  10  hours  and  30  minutes,  but  the  length  of  time  the  aircraft  was  actually 
airborne  ranged  from  30  minutes  to  nine  hours  and  30  minutes.  This  research  only 
evaluated  Global  Hawk  data  during  its  flight,  so  information  provided  when  it  was  on  the 
ground  was  disregarded. 


-118.6-118.2  -117.8  -117.4  -117-116.7 
PRIMARY  WGS84  LONG 


Figure  1 1 .  Aeronautical  map  of  Edwards  AFB  military  operating  areas.  The  red  dot 
indicates  the  Edwards  AFB  rawinsonde  site  (Maptech  2003). 
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PRIMARY  WGS84  LONG 


Figure  12.  Latitude  versus  longitude  plot  of  Global  Hawk  flight  routes  for  all  five 
examined  missions.  The  red  dot  indieates  the  Edwards  AFB  rawinsonde  site.  The 
latitude/longitude  range  eorresponds  to  Figure  11. 


The  Air  Data  System  onboard  the  Global  Hawk  compiles,  processes,  and 
transmits  all  atmospheric  variables  to  the  Integrated  Mission  Management  Computers 
(IMMC)  onboard  the  aircraft.  These  parameters  include  static  air  pressure,  outside  air 
temperature,  wind  speed,  and  wind  direction.  The  units  of  measurement  for  each  of  these 
variables  are  pounds  per  square  foot,  degrees  Celsius,  degrees  true,  and  knots; 
respectively  (Northrop  Grumman  2001). 


2)  ASPAM 

AFCCC  manages  an  extensive  database  of  environmental  information  which  is 
used  as  input  for  ASPAM.  The  input  consists  of  surface  observations,  upper-air 
soundings,  satellite  observations,  snow  and  ice  field  data,  solar  and  geomagnetic 
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information,  aerosol  information,  Digital  Terrain  Elevation  Data  (DEED),  Real-Time 
Nephanalysis  (RTNEPH)  data,  NOGAPS  information,  aireraft  reports  (AIREPS),  pilot 
reports  (PIREPS),  and  Aireraft  Communieation  Addressing  and  Reporting  System 
(ACARS)  reports.  The  profiles  provide  data  from  the  surfaee  to  100,000  feet  in  500-foot 
inerements  for  specified  latitudes,  longitudes,  and  times.  Among  the  data  in  the  profiles 
are  pressure,  temperature,  wind  direction,  and  wind  speed  in  units  of  millibars,  degrees 
Celsius,  degrees  true,  and  meters  per  second;  respectively.  ASPAM  uses  the  method  of 
Multivariate  Optimal  Interpolation  (MVOI)  to  determine  data  values  at  each  level. 

MVOI  assigns  a  value  of  significance  to  each  input  data  point  based  on  its  distance  from 
the  profile  location,  the  difference  between  the  observation  time  and  the  profile  time,  the 
accuracy  of  the  observing  equipment,  and  the  expected  accuracy  of  the  first  guess.  Then 
all  of  the  weighted  inputs  are  combined  to  create  a  “best  guess”  vertical  profile  at  the 
given  position  and  time.  Two  ASPAM  vertical  profiles  per  Global  Hawk  mission  were 
obtained  for  this  research,  one  for  the  ascending  segment  and  the  other  for  the  descending 
segment  of  the  flight  profiles. 

3)  RAOB 

The  weather  unit  at  Edwards  AEB  supports  the  flight  testing  squadrons  by 
providing  surface  and  upper-air  observations,  as  well  as  forecasts.  The  unit  launches 
rawinsondes  based  on  customer  requirements  and  the  corresponding  upper-air  data  are 
posted  to  its  web  site  and  archived  for  future  use.  Data  were  available  for  three  of  the 
five  Global  Hawk  mission  dates  within  four  hours  of  the  flight  times.  Another  dataset 
was  available  for  one  of  the  other  mission  dates  but  was  approximately  15  to  18  hours 


22 


prior  to  the  flight  time.  This  dataset  was  not  used  beeause  it  was  not  eonsidered  to  be 
representative  of  the  atmospherie  eonditions  during  the  Global  Hawk  flight.  These 
upper-air  datasets  were  obtained  in  text  format  and  eompared  with  the  Global  Hawk 
mission  data  as  a  seeond  test  of  aeeuraey. 

b.  Data  limitations 

Errors  and  limitations  are  inherent  in  any  atmospherie  proeess  or  dataset.  The 
results  of  this  researeh  are  of  value,  so  long  as  these  problems  are  identified  and  taken 
into  aeeount. 

1)  Global  Hawk 

Data  reeeived  from  the  Global  Hawk  Air  Data  System  (ADS)  are  oeeasionally 
impreeise  or  ineorreetly  formatted  during  transmission.  The  ADS  is  a  redundant  system, 
whieh  means  another  ADS  proeessor  relays  data  to  a  separate  Integrated  Mission 
Management  Computer  (IMMC).  The  Global  Hawk  is  also  equipped  with  a  system 
status  eomputer,  whieh  reports  any  expeeted  problems  in  ADS  output.  Distorted 
information  ean  still  make  it  through  the  system,  however,  and  all  datasets  must  be 
eheeked  for  erroneous  or  missing  data.  This  researeh  will  help  determine  if  and  when 
ADS  output  eannot  be  trusted  as  an  aeeurate  data  souree. 

2)  ASPAM 

ASP  AM  interpolates  data  to  a  speeifie  loeation  and  time  and  the  profiles  are  a 
representation  of  the  atmosphere  based  on  the  input.  They  may  or  may  not  refieet  the 
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actual  conditions.  Therefore,  any  errors  in  the  input  will  result  in  less  accurate  ASPAM 
profiles.  Possible  eauses  for  input  error  include  instrumentation  faults,  problems  in  data 
entry  or  transmission,  and  inconsistencies  in  equipment  and  reporting. 

3)  RAOB 

Possible  errors  in  rawinsonde  data  include  transmission  problems,  incorreet  data 
entry,  and  equipment  faults.  The  upper- air  datasets  used  in  this  researeh  did  not  exactly 
coincide  with  the  Global  Hawk  mission  times  or  locations,  but  were  as  close  as  feasibly 
possible  to  perform  statistieal  analyses  and  comparisons.  The  results  of  these  tests  can  be 
used  to  support  the  evaluations  based  on  ASPAM  data  and  gain  further  insight  into  the 
validity  of  the  Global  Hawk  data.  Allen  (2003)  states,  “the  raw  sounding  data  includes 
the  temperature,  dew  point  temperature,  pressure,  and  wind  speed  and  direction  for  the 
atmospheric  column.  This  column  is  assumed  to  be  directly  over  the  station  that  launehes 
the  balloon.  However,  this  assumption  is  not  entirely  accurate  since  the  atmospheric  flow 
will  typically  carry  the  balloon  downwind  away  from  the  station  as  it  ascends.  Errors  in 
balloon  location  and  height  can  occur  if  the  balloon  is  not  rising  at  the  standard  ascension 
rate  of  300  meters  per  minute.  Errors  in  sounding  data  ean  also  oecur  due  to  the  lag  time 
of  sensors  in  the  instrumentation  paekage.  Eor  example,  if  the  balloon  is  passing  through 
a  rapidly  ehanging  layer,  some  of  the  information  will  not  get  reported  since  the 
instruments  cannot  react  to  the  changes  fast  enough.  The  typical  lag  for  temperature 
sensors  is  4  to  20  seconds,  with  the  lag  increasing  with  altitude.  The  pressure 
measurements  used  for  determining  the  height  have  a  standard  error  of  ±1  mb  at  the 
surface  and  10  mb  with  an  error  of  2  mb  at  500  mb  (Golden  et  al.  1986).” 
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Edwards  AFB  is  located  in  the  Mojave  Desert  and  the  weather  in  this  region  is 
generally  stable,  espeeially  during  the  late  spring  and  summer.  Upper-air  eharts  were 
obtained  eoineident  to  the  five  Global  Hawk  missions  and  it  was  determined  that  the 
weather  eonditions  were  stable  during  these  times.  Upper-level  wind  speeds  during  the 
first  three  missions  were  such  that  the  rawinsondes  may  not  have  drifted  signifieantly  far 
downwind.  Considerably  stronger  upper-level  wind  speeds  and  substantial  direetional 
shear  during  the  two  August  missions  may  have  caused  the  rawinsondes  to  drift  farther 
downstream  produeing  RAOBs  that  were  not  representative  of  the  atmospheric 
conditions  in  the  Global  Hawk  flight  areas.  See  Appendix  B  for  the  upper-level  wind 
analyses  corresponding  to  all  five  mission  dates. 
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4,  Methodology 

a.  Data  acquisition 

Three  sets  of  data  were  analyzed  and  eompared  in  this  researeh.  Environmental 
information  from  the  Global  Hawk  was  the  foeus  of  this  researeh  and  was,  therefore,  the 
primary  dataset  examined.  The  dates,  times,  and  loeations  of  the  Global  Hawk  data 
determined  the  format  of  the  additional  data  requested.  The  additional  data  ineluded 
vertical  atmospheric  profiles  from  the  Atmospheric  Slant  Path  Analysis  Model  (ASPAM) 
and  upper-air  observations  from  Edwards  AFB  rawinsondes. 

1)  Global  Hawk 

Global  Hawk  environmental  information  was  obtained  from  flight-test  databases 
maintained  by  the  452”‘*  FLTS  at  Edwards  AFB.  The  available  weather  data  included 
pressure,  temperature,  wind  direction,  and  wind  speed  in  units  of  pounds  per  square  foot, 
degrees  Celsius,  degrees  true,  and  knots;  respectively.  Aircraft  flight  parameters  were 
collected  along  with  the  weather  information,  including  altitude  in  feet  above  mean  sea 
level  (MSF),  ascent  mode,  descent  mode,  altitude  hold  mode,  pitch  angle,  and  roll  angle. 
This  information  was  used  to  determine  if  the  Global  Hawk  environmental  data  were 
influenced  by  various  aspects  of  the  aircraft’s  flight.  The  data  reports  were  available 
once  per  second  for  each  mission  but  only  data  at  specific  flight  levels  were  analyzed  and 
compared.  Data  from  five  separate  test  missions  of  the  same  aircraft.  Air  Vehicle  7 
(AV7),  were  provided  in  comma-separated  format  on  CD-ROM  media. 
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2)  ASPAM 

AFCCC  provided  vertical  atmospheric  profiles  from  the  Atmospheric  Slant  Path 
Analysis  Model  (ASPAM).  The  profiles  provided  weather  information  from  the  surface 
to  100,000  feet  MSL  in  500-foot  increments.  The  environmental  data  used  in  this 
research  included  height,  pressure,  temperature,  wind  direction,  and  wind  speed,  in  units 
of  feet  MSL,  millibars,  degrees  Celsius,  degrees  true,  and  meters  per  second; 
respectively. 

A  total  of  ten  ASPAM  profiles  were  requested  from  AFCCC,  one  for  the 
ascending  and  one  for  the  descending  portion  of  each  Global  Hawk  mission.  The 
parameters  of  the  ASPAM  profiles  were  based  on  the  median  times,  latitudes,  and 
longitudes  for  each  ascent  and  descent  of  the  Global  Hawk  missions.  The  profiles  were 
then  created  and  e-mailed  in  text  format. 

3)  RAOB 

Rawinsonde  observations  were  provided  by  the  Edwards  AFB  weather  unit.  The 
observations  provided  altitude,  pressure,  temperature,  wind  direction,  and  wind  speed  in 
units  of  feet  MSL,  millibars,  degrees  true,  and  knots;  respectively.  RAOBs  were 
available  on  three  of  the  Global  Hawk  mission  dates  within  four  hours  of  the  flight  times. 
An  additional  RAOB  was  available  for  another  mission  but  was  15  to  18  hours  prior  to 
the  Global  Hawk  flight  time  and  was,  therefore,  not  used  in  this  research.  All 
Rawinsonde  observations  were  compiled  and  e-mailed  in  text  format. 
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b.  Data  formatting 

Global  Hawk,  ASP  AM,  and  RAOB  files  were  imported  into  Microsoft  Excel  and 
SAS  Institute  JMP  and  saved  in  the  same  file  format  prior  to  quality  checks.  All  of  the 
datasets  were  then  checked  for  errors  and  any  faulty  reports  were  removed  prior  to 
statistical  analysis.  All  three  data  types  were  arranged  in  the  same  format  and  unit 
convention  to  simplify  comparisons  and  tests. 

1)  Global  Hawk 

Comma-separated  text  files  were  formatted  and  imported  into  statistical  software. 
Datasets  from  all  five  Global  Hawk  missions  were  carefully  examined  for  obvious  errors 
or  incorrect  entries.  If  any  variable  in  a  row  of  data  was  invalid  or  suspect,  the  entire  row 
of  data  was  removed.  Additionally,  all  data  reports  while  the  aircraft  was  on  the  ground 
were  excluded  from  analysis  because  this  research  focused  only  on  verifying  data  during 
ascent  and  descent.  Global  Hawk  data  are  reported  once  per  second  but  the  only  data  of 
interest  were  at  altitudes  from  the  surface  (approximately  2,500  feet  MSL)  to  maximum 
flight  level  (20,000  to  60,000  feet  MSL)  in  500-foot  increments.  Therefore,  all  data 
reports  were  removed  except  for  those  closest  to  a  specified  500-foot  flight  level.  Table  2 
depicts  a  sample  of  the  Global  Hawk  data  analyzed  and  the  difference  between  the  actual 
altitude  and  the  corresponding  500-foot  level.  The  table  data  represent  the  greatest 
differences  for  the  ascending  and  descending  portions  of  each  mission.  Appendix  C 
provides  a  complete  listing  of  the  entire  Global  Hawk  dataset  which  was  used  in  this 
research. 
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Table  2.  Global  Hawk  flight  levels  examined  and  differences  from  corresponding  500- 
foot  levels.  Data  displayed  below  represent  the  levels  with  the  greatest  differences  per 
mission.  See  Appendix  C  for  the  complete  dataset. 


Global  Hawk 
Mission  # 

Global  Hawk 
Flight  Mode 

Global  Hawk 
Flight  Level 
(ft  MSL) 

500-Foot 
Flight  Levels 
(ft  MSL) 

Difference 
(ft  MSL) 

1 

Ascent 

15972.29 

16000 

27.71 

1 

Descent 

18528.68 

18500 

28.68 

2 

Ascent 

5028.65 

5000 

28.65 

2 

Descent 

60002.20 

60000 

2.20 

3 

Ascent 

7527.80 

7500 

27.80 

3 

Descent 

27969.59 

28000 

30.41 

4 

Ascent 

12028.23 

12000 

28.23 

4 

Descent 

50479.88 

50500 

20.12 

5 

Ascent 

14028.72 

14000 

28.72 

5 

Descent 

55038.77 

55000 

38.77 

The  Global  Hawk  continuously  checks  the  Air  Data  System  (ADS)  and  indicates 
any  suspected  problems  in  the  ADS  Status  column.  Thusly,  any  reports  with  ADS 
problems  were  removed  from  the  datasets  prior  to  analysis.  A  total  of  772  rows  of 
Global  Hawk  data  were  available  for  analysis  once  all  invalid  or  suspect  rows  were 
removed. 

Global  Hawk  reports  of  Climb  Mode,  Altitude  Hold  Mode,  and  Descend  Mode 
were  checked  for  errors  and  some  possible  system  faults  were  identified.  For  example, 
the  aircraft  continued  to  climb,  sometimes  as  much  as  2,500  feet,  while  still  reporting  to 
be  in  Altitude  Hold  Mode.  During  other  missions,  the  aircraft  indicated  to  be  in  Altitude 
Hold  Mode  but  its  altitude  fluctuated  as  much  as  6,000  feet.  Most  of  these  suspect  data 
reports  were  not  included  in  the  analysis  but  some  of  them  had  to  be  kept  to  ensure  a 
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complete  mission  profile  was  examined.  Notes  pertaining  to  the  Global  Hawk  data  and 
varying  flight  modes  are  shown  in  Appendix  D. 

The  wind  data  reported  by  the  Global  Hawk  were  not  in  the  standard 
meteorologieal  format  and  had  to  be  eonverted  prior  to  analysis.  Original  values  of  wind 
direetion  were  between  -180  and  180  degrees  and  depleted  winds  toward,  rather  than 
from  the  direetion.  It  was  also  determined  that  negative  values  eorresponded  to  winds  in 
the  western  quadrants  (west  of  north)  and  positive  values  eorresponded  to  winds  in  the 
eastern  quadrants  (east  of  north).  The  method  used  to  put  the  Global  Hawk  wind 
direetion  into  the  standard  meteorologieal  format  was  to  add  180  degrees  to  all  values. 
The  wind  speed  data  were  reported  in  units  of  knots  and  were  eonverted  to  meters  per 
seeond  by  dividing  all  values  by  1.943. 

The  pressure  data  from  the  Global  Hawk  were  eonverted  from  pounds  per  square 
foot  to  millibars  by  dividing  all  values  by  2.08854.  This  ensured  all  pressure  values  were 
eonsistently  formatted  prior  to  eomparison. 

2)  ASPAM 

All  ASPAM  data  files  were  eonverted  into  Exeel  and  JMP  formats  and  eombined 
with  the  Global  Hawk  data.  The  vertieal  profiles  were  ereated  in  the  standard 
meteorologieal  format;  therefore,  no  unit  eonversions  were  neeessary.  All  applieable 
data  were  reported  from  the  surfaee  to  100,000  feet  MSL  in  500-foot  inerements  but  only 
data  with  eorresponding  Global  Hawk  reports  were  examined.  Therefore,  772  rows  of 
ASPAM  data  were  available  for  analysis  and  eomparison  with  the  Global  Hawk  data. 
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3)  RAOB 

All  Edwards  AFB  rawinsonde  observations  were  converted  into  Excel  and  JMP 
formats  and  combined  with  Global  Hawk  and  ASPAM  data.  Any  additional  information, 
such  as  density  and  humidity,  was  removed  before  incorporating  the  observations  with 
the  other  datasets.  Wind  speed  was  reported  in  knots  and  was  converted  to  meters  per 
second  by  dividing  all  values  by  1 .943.  Only  three  RAOBs  were  used  in  this  research, 
the  dates  and  times  of  which  were  29  May  03  at  163 IZ,  8  Aug  03  at  1230Z,  and  15  Aug 
03  at  1900Z.  RAOB  data  were  provided  in  500-foot  increments  from  the  surface  to 
15,000  feet  and  then  in  1,000-foot  increments  from  15,000  feet  to  the  end  of  the 
transmission.  The  final  altitudes  reported  by  each  RAOB  were:  29  May  -  40,000  feet,  8 
Aug  -  57,000  feet,  and  15  Aug  -  48,000  feet.  Available  RAOB  data  points  were 
significantly  less  than  Global  Hawk  and  ASPAM  reports  due  to  missing  RAOBs  for 
missions  two  and  three,  fewer  RAOB  reports  above  15,000  feet,  and  lower  maximum 
altitudes  of  the  balloon  datasets.  Therefore,  246  rows  of  RAOB  data  were  analyzed  and 
compared  with  the  Global  Hawk  and  ASPAM  data. 

c.  Data  verification 

Once  all  of  the  Global  Hawk,  ASPAM,  and  RAOB  datasets  were  formatted  and 
standardized,  several  statistical  tests  were  accomplished  to  verify  the  Global  Hawk  data. 
The  primary  test  performed  was  simple  linear  regression  of  Global  Hawk  data  with 
ASPAM  and  RAOB  data. 
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1)  Statistical  OVERVIEW 

Regression  analysis,  as  stated  by  Montgomery  and  Runger  (2003),  is  a  useful 
method  to  determine  the  relationship  between  two  or  more  variables.  Simple  linear 
regression  was  used  to  resolve  relationships  between  the  Global  Hawk,  ASPAM,  and 
RAOB  variables  of  pressure,  temperature,  wind  speed,  and  wind  direetion.  Global  Hawk 
data  were  used  as  the  independent,  or  predictor  variables  represented  by  x  and  ASPAM 
and  RAOB  data  were  used  as  the  dependent  variables,  or  predictands  represented  by  y 
(Wilks  1995).  Regression  analysis  provided  valuable  information  on  the  accuracy  of  the 
Global  Hawk  data  after  plotting  corresponding  datasets  on  a  scatter  diagram.  For 
example,  a  plot  was  created  of  Global  Hawk  pressure  (x)  versus  ASPAM  pressure  (y) 
(Figure  13)  where  each  (x,  ,  y/)  pair  was  shown  as  a  single  point  on  the  graph.  The  data 
points  appeared  to  be  randomly  scattered  around  a  straight  line  and  a  simple  linear 
regression  model  was  developed.  The  model  is  given  by  the  equation 

Y  =  (1) 

where  (3^  and  f3^  are  unknown  regression  coefficients  representing  the  intercept  and 
slope  of  the  linear  model,  respectively.  The  random  errors,  e  ,  were  assumed  to  be 
uncorrelated  random  variables  with  a  mean  of  zero  and  unknown  variance  <j^ 
(Montgomery  and  Runger  2003). 
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GH  Static  Pressure  (mb) 


Figure  13.  Scatterplot  of  Global  Hawk  pressure  (x)  versus  ASPAM  pressure  (y). 


This  model  was  used  to  determine  relationships  between  Global  Hawk  variables  and 
ASPAM  and  RAOB  variables.  In  the  words  of  an  Air  Force  Institute  of  Technology 
(AFIT)  statistics  professor,  “There  are  good  models  and  bad  models,  but  no  correct  ones. 
It’s  a  model.  If  it  were  right,  it  would  be  scientific  fact  (R.N.  Benton  2003,  personal 
communication).”  The  estimates  of  and  produce  a  “best  fit”  line  to  the  data  and 

the  technique  of  least  squares  regression  is  performed  to  minimize  the  sum  of  the  squares 
of  the  vertical  deviations,  or  total  spread,  of  the  y-values  from  the  line.  The  least  squares 
estimates  of  the  intercept  and  slope  in  the  linear  regression  model  with  n  observations  are 

k  =  y-k^  (2) 

and 


X:  - 


f 

V  1=1 


V  i=i  y 


V  1=1  y 


n 


(3) 
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where 


-  Ex 

\n)tt 


(4) 


and 


X  = 


\n) 


(5) 


The  numerator  and  denominator  of  Equation  3  are  abbreviated  as  and  SS^^, 


respectively.  The  expression  SS^  measures  the  spread  of  x„  which  is  the  sum  of  squares 
around  the  mean,  x  .  Another  useful  term  is  SS  ,  which  measures  the  spread  around  y  , 


and  is  given  by 


(6) 


/=1 


This  value  is  also  identified  as  the  total  sum  of  squares,  SSj, ,  and  is  used  to  test  the 
significance  of  regression  for  each  model.  The  regression  line  or  “best  fit”  becomes 

y  =  (7) 

Each  {xi,  yt)  pair  of  observations  fit  the  relationship 

(8) 


where  is  called  the  error  or  residual,  which  indicates  the  vertical  distances  between  the 
points  and  the  fitted  line.  This  term  is  given  by 

(9) 

The  residuals  provide  information  on  the  adequacy  of  each  linear  model  and  are  used  to 
determine  an  estimate  of  the  corresponding  variance,  <j^ .  The  analysis  of  variance 
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(ANOVA)  table  was  used  to  determine  the  value  of  eaeh  regression  model.  The  analysis 
of  varianee  identity  is 

(10) 

i=\  i=\  i=\ 

The  first  value  on  the  right  side  of  the  equation  measures  the  amount  of  variability  iny, 
aeeounted  for  by  the  regression  line  and  is  referred  to  as  the  regression  sum  of  squares, 
S'5'^ .  The  seeond  term  on  the  right  side  of  the  equation  measures  the  residual  variation 

left  unexplained  by  the  regression  line  and  is  referred  to  as  the  error  sum  of  squares,  SS^ . 
It  is  used  to  show  how  mueh  the  predieted  values,  y . ,  differ  from  the  true  y-values  and  is 
defined  as 

S'S£=Z(.'’,-v,f  (11) 

i=\  i=\ 

Several  statisties  provided  information  on  the  “goodness  of  fit”  of  the  models. 

The  first  statistie  of  interest  is  the  mean  square  error,  MSB,  whieh  is  an  estimator  of  cr^ 
and  represents  the  sample  varianee  of  the  residuals.  It  is  ealeulated  by  dividing  the  error 
sum  of  squares  by  its  number  of  degrees  of  freedom  (Montgomery  and  Runger  2003), 
whieh  is  shown  by 

MSE  =  (j^  (12) 

n-2 
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The  root  mean  square  error,  RMSE,  is  the  square  root  of  the  mean  square  error 
and  is  used  beeause  it  retains  the  units  of  the  foreeast  variable.  The  magnitude  of  the 
error  is,  therefore,  easier  to  interpret  using  RMSE  (Wilks  1995). 

A  second  statistic  used  to  judge  the  adequacy  of  the  regression  models  is  the 
coefficient  of  determination,  R^,  which  is  defined  by 


SSj  SSj  ■ 


(13) 


It  is  the  proportion  of  the  variation  of  the  predictand  that  is  “accounted  for”  by  the 
regression.  A  perfect  linear  regression  model  would  have  =  \  because  SSj^  =  SSj  and 


SSj,  =  0  and  a  completely  useless  model  would  have  =  0  because  =  0  and 
SSj,  =  SSj, .  Thus,  higher  values  of  R^  indicate  a  more  valuable  model  (Wilks  1995). 

Another  useful  statistic  is  the  F-ratio,  which  is  listed  on  the  ANOVA  table  and  is 
used  to  describe  the  strength  of  regression  models.  It  is  the  ratio  of  the  regression  sum  of 
squares,  SSj^ ,  with  the  mean  square  error,MSE.  Higher  values  of  F  indicate  more 

effective  regression  because  a  strong  relationship  between  x  and  y  produces  a  large  SSj^ 
and  a  small  MSE  (Wilks  1995). 

Once  simple  linear  regression  was  accomplished  between  each  Global  Hawk 
variable  and  the  corresponding  ASPAM  and  RAOB  variables,  hypothesis  tests  were 
performed  to  further  justify  each  model.  These  tests  were  performed  to  determine  if  the 
slope  of  each  regression  line  was  equal  to  one,  =  1 .  If  the  hypothesis  tests  passed,  then 
the  relationships  between  the  Global  Hawk  variables,  x,  and  the  ASPAM  and  RAOB 
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variables,  y,  could  be  validated.  The  Mest  was  employed  to  test  the  hypothesis  of  =  1 


for  each  regression  model.  The  test-statistic  defined  as 


A-a,o 


(14) 


follows  the  t-distribution  with  n-2  degrees  of  freedom  where  o  =  1  for  the  tests  used 


in  this  research.  The  hypotheses  of  o  =  1  is  rejected  if  [tol  >  n-i  ’  where 

a  =  0.05  and  n  represents  the  number  of  observations  used  in  each  model.  Rejection  of 
the  null  hypothesis  occurs  if  the  probability  of  observing  the  attained  or  more  extreme 
value  of  the  test  statistic,  given  the  null  hypothesis  is  true,  or  p-value,  is  less  than  or  equal 
to  the  test  level  of  «  =  0.05 .  If  the  hypothesis  tests  were  rejected  then  it  could  not  be 
concluded  that  the  slopes  of  the  regression  lines  were  equivalent  to  one.  The  hypothesis 
tests  of  y^i  =  1  used  in  this  study  were  likely  to  fail  because  each  of  the  linear  models  in 
this  research  had  an  extremely  small  variance.  Models  with  low  variance  tend  to  have 
large  values  of  to  based  on  the  definition  of  the  test  statistic  in  Equation  14.  The 
regression  lines  in  this  research  were  not  invalidated  based  on  the  results  of  these 
hypothesis  tests,  however,  because  the  variances  were  low  and  the  fit  was  good  for  all 
models. 

Another  hypothesis  test  used  to  illustrate  the  significance  of  regression  was 
performed  by  testing  if  the  slopes  of  the  regression  lines  were  zero,  y^i  =  0 .  If  a  test 

failed,  meaning  the  hypothesis  was  rejected,  then  it  was  concluded  that  there  was  a 
significant  relationship  between  x  andy.  Rejection  of  the  null  hypothesis  occurs  if  the 
probability  of  observing  the  attained  or  more  extreme  value  of  the  test  statistic,  given  the 
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null  hypothesis  is  true,  or  p-value,  is  less  than  or  equal  to  the  test  level  of  a  =  0.05  . 
Rejecting  the  hypothesis  implies  that  x  is  of  value  in  explaining  the  variability  of  y  and 
could  mean  that  the  linear  model  is  adequate.  Failure  to  reject  the  hypothesis,  =  0  , 

means  thatx  is  of  little  value  in  explaining  the  variability  iny  and  is  the  same  as 
concluding  that  there  is  no  linear  relationship  between  x  andy  (Montgomery  and  Runger 
2003).  The  tests  of  =  0  used  in  this  research  concluded  0 .  Failure  of  these  tests 
implied  that  the  Global  Hawk  variables  were  of  use  in  determining  the  “ground  truth” 
variables  of  ASPAM  and  RAOB  and  that  there  was  a  linear  relationship  between  them. 


2)  Wind  direction  techniques 

Wind  direction  verification  was  difficult  due  to  the  non-linearity  of  the  data.  A 
method  of  validating  wind  data  used  in  meteorology  is  the  root-mean-square  or  RMS 
vector  error.  The  wind  data  were  first  broken  down  into  their  respective  easterly  (u)  and 
northerly  (v)  components  by 

u  =  -Xsm0  (15) 

and 


v  =  -Xcos^  (16) 

where  X  is  the  magnitude  of  the  wind  speed  in  m  s'^  and  0  is  the  direction  of  the  wind  in 
degrees  measured  clockwise  from  north.  The  wind  components  were  obtained  for  all 
three  datasets  and  the  RMS  vector  errors  were  calculated  for  Global  Hawk  versus 
ASPAM  and  Global  Hawk  versus  RAOB.  The  RMS  vector  error  is  defined  as 


n 


(17) 
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where  V  represents  the  magnitude  of  the  RMS  veetor  error,  n  is  the  total  number  of 
observations,  u i  is  the  Global  Hawk  easterly  eomponent,  U2  is  the  ASPAM  or  RAOB 
easterly  eomponent,  v/  is  the  Global  Hawk  northerly  eomponent,  and  V2  is  the  ASPAM  or 
RAOB  northerly  eomponent. 

Another  method  used  to  verify  wind  direetion  was  to  determine  the  mean 
differenees  between  Global  Hawk  and  “ground  truth”  data.  Two  eolumns  were  ereated  in 
JMP  to  ealeulate  the  differenees  of  Global  Hawk  wind  direetions  with  ASPAM  and 
RAOB  wind  direetions.  A  formula  was  used  for  eaeh  column  to  guarantee  the  acute 
angle  was  given  for  the  difference  between  the  wind  directions.  The  formula  used  in 
JMP  returned  the  minimum  of  three  calculations,  which  were 


(18) 

A1-A2+360 

(19) 

A1-A2-360 

(20) 

where  Xj  is  the  Global  Hawk  direction  and  X2  is  the  ASPAM  or  RAOB  direction. 

Wind  direction  information  was  also  examined  as  circular  statistical  data.  The 
distributions  of  wind  direction  reports  could  not  be  considered  normal  because  the  data 
was  non-linear.  The  von  Mises  distribution  for  circular  data  is  analogous  to  the  normal 
distribution  for  linear  data  (Fisher  1993).  A  statistical  program  called  Oriana  was  used  to 
analyze  the  circular  datasets  and  the  wind  direction  reports  were  shown  to  fit  the  von 
Mises  distribution.  Additional  analyses  using  circular  statistical  information  were 
possible  but  time  constraints  did  not  allow  further  investigation  into  this  area. 
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5,  Analysis  and  Results 

a.  Pressure 

All  tests  and  analyses  of  Global  Hawk  pressure  were  very  successful  and  provided 
sufficient  information  to  conclude  that  Global  Hawk  pressure  reports  were  accurate.  The 
linear  regression  of  Global  Hawk  pressure  with  ASPAM  and  RAOB  pressure  produced 
significant  linear  models  and  extremely  low  values  of  variance  and  RMSE. 


1)  ASPAM 

All  five  Global  Hawk  missions  provided  772  data  points  used  in  linear  regression 
of  Global  Hawk  pressure  (x)  with  ASPAM  pressure  (y).  The  IMP  scatterplot  and  fitted 
line  is  depicted  in  Figure  14  and  the  regression  analysis  results  are  illustrated  in  Figure 
15.  The  results  include  R^,  RMSE,  ANOVA  table,  F-ratio,  and  parameter  estimates, 
along  with  the  t-test  values  for  the  hypothesis  test  =  0  .  The  figures  show  the  Global 
Hawk  is  an  accurate  predictor  of  atmospheric  pressure  as  established  by  ASPAM. 


ASPAM  Pressure  (mb)  Predicted  PO.OOOO 
RSq=1.00  RMSE=8.2911 


Figure  14.  Simple  linear  regression  scatterplot  and  fitted  line  of  Global  Hawk  pressure 
(x)  versus  ASPAM  pressure  (y). 
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Summary  of  Fit 

RSquare 

0.99867 

RSquare  Adj 

0.998668 

Root  Mean  Square  Error 

8.29109 

Mean  of  Response 

347.1323 

Observations  (or  Sum  Wgts) 

772 

Analysis  of  Variance 

F  Ratio 
578249.4 
Prob  >  F 
0.0000 


Source  DF  Sum  of  Squares  Mean  Square 

Model  1  39750116  39750116 

Error  770  52931  68.742169 

C.  Total  771  39803048 


Parameter  Estimates 

Term  Estimate  Std  Error  t  Ratio  Prob>|t| 

Intercept  9.1112201  0.535386  17.02  <.0001 

GH  Static  Pressure  (mb)  0.9157469  0.001204  760.43  0.0000 

Figure  15.  Regression  analysis  results  of  Global  Hawk  pressure  (x)  versus  ASPAM 
pressure  (y). 


The  hypothesis  test  for  \  was  accomplished  using  the  JMP  output  from 
Figure  15  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
freedom  are  defined  by  n  -  2  =  770  for  Global  Hawk  versus  ASPAM  regression.  The 
JMP  values  used  were  =  0.9157469  and  =  68.742169 .  The  critical  t-value  and 
were  calculated  as  =^0025770  =1-963  and  =47400996.86.  The  output 
t-value  was  =  -69.963  which  led  to  the  conclusion  of  rejecting  the  null  hypothesis 
=  1 .  As  previously  stated,  this  test  was  likely  to  fail  because  of  the  extremely  low 
variance  of  the  linear  model.  Therefore,  all  assumptions  of  Global  Hawk  pressure  data 
quality  were  based  on  the  regression  analysis  and  hypothesis  test  of  =  0  . 
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Residual  analysis  showed  there  were  no  significant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  scale  as  the  regression  variables. 
Close  examination  of  the  data  revealed  some  possible  patterns  in  the  bias,  however. 
Global  Hawk  pressure  values  primarily  had  a  positive  bias  during  ascent  and  a  negative 
bias  during  descent.  Also,  the  amount  of  bias  seemed  to  change  in  relation  to  the  rate  at 
which  the  Global  Hawk  climbed  or  descended.  It  is  possible  that  a  time  delay  between 
measuring  and  transmitting  atmospheric  values  by  the  Global  Hawk  would  induce  the 
bias  observed  in  this  study.  During  ascent  (descent),  the  Global  Hawk  could  measure  the 
pressure  at  an  altitude  lower  (higher)  than  the  height  at  which  it  transmitted  the  data. 

This  slight  difference  in  flight  level  could  lead  to  pressure  values  higher  (lower)  than 
those  at  which  the  data  was  relayed.  This  height  difference  could  vary  in  relation  to  the 
rate  at  which  the  Global  Hawk  climbed  or  descended  and  the  pressure  bias  would  reflect 
this  difference.  A  possible  reason  for  this  slight  difference  could  be  the  vertical 
separation  between  the  weather  sensors  at  the  rear  of  the  aircraft  and  the  altimeter  during 
ascent  and  descent.  See  Appendix  E  for  an  illustration  of  the  Global  Hawk  dimensions 
and  location  of  the  sensors.  The  residual  versus  predicted  plots  for  each  Global  Hawk 
mission  are  depicted  in  Figure  16. 
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ASPAM  Pressure  (mb)  Residual  ASPAM  Pressure  (mb)  Residual 


(a) 


(b) 


ASPAM  Pressure  (mb)  Predicted 


(c) 


ASPAM  Pressure  (mb)  Predicted 


(d) 


ASPAM  Pressure  (mb)  Predicted 


(e) 


ASPAM  Pressure  (mb)  Predicted 


Figure  16.  Residual  versus  predicted  plots  of  ASPAM  pressure  for  Global  Flawk 
missions  (a-e)  1  through  5. 
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2)  RAOB 


Only  three  rawinsonde  observations  were  available  for  comparison  with  Global 
Hawk  data  and  a  total  of  246  reports  were  analyzed.  Global  Hawk  pressure  was  used  as  a 
predictor  for  RAOB  pressure  and  the  scatterplot  with  fitted  line  is  shown  in  Figure  17. 

■y 

The  JMP  regression  analysis  results  are  shown  in  Figure  18,  including  R  ,  RMSE, 
ANOVA  table,  F-ratio,  and  parameter  estimates  including  the  t-test  results  for  the 
hypothesis  =  0  .  Regression  analysis  showed  the  Global  Hawk  is  an  accurate  predictor 
of  atmospheric  pressure  as  established  by  RAOBs. 


100  200  300  400  500  600  700  800  900 
RAOB  Pressure  (mb)  Predicted  PO.OOOO 
RSq=1.00  RMSE=5.5163 


Figure  17.  Simple  linear  regression  scatterplot  and  fitted  line  of  Global  Hawk  pressure 
(x)  versus  RAOB  pressure  (y). 

The  hypothesis  test  for  \  was  accomplished  using  the  JMP  output  from 

Figure  18  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
freedom  are  defined  by  n-2  =  244  for  Global  Hawk  versus  RAOB  regression.  The  JMP 
values  used  were  =  0.9874827  and  =  30.429904  .  The  critical  t-value  and  SS^ 
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were  calculated  as  =^0025244  =1-9697  and  SS^  =47400996.86.  The  output  t- 
value  was  =  -15.263  which  led  to  the  conclusion  of  rejecting  the  null  hypothesis 
=  1 .  As  previously  stated,  this  test  was  likely  to  fail  because  of  the  extremely  low 
variance  of  the  linear  model.  Therefore,  all  assumptions  of  Global  Hawk  pressure  data 
quality  were  based  on  the  regression  analysis  and  hypothesis  test  of  =  0  . 


Summary  of  Fit 


RSquare 

0.999509 

RSquare  Adj 

0.999507 

Root  Mean  Square  Error 

5.516331 

Mean  of  Response 

517.4063 

Observations  (or  Sum  Wgts) 

246 

Analysis  of  Variance 

F  Ratio 
496307.6 
Prob  >  F 
0.0000 


Source  DF  Sum  of  Squares  Mean  Square 

Model  1  15102592  15102592 

Error  244  7425  30.429904 

C.  Total  245  15110017 


Parameter  Estimates 

Term  Estimate  Std  Error  t  Ratio  Prob>|t| 

Intercept  23.2148  0.784719  29.58  <.0001 

GH  Static  Pressure  (mb)  0.9874827  0.001402  704.49  0.0000 

Figure  18.  Regression  analysis  results  of  Global  Hawk  pressure  (x)  versus  RAOB 
pressure  (y). 

Residual  analysis  showed  there  were  no  significant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  scale  as  the  regression  variables. 
Similar  patterns  in  the  bias  were  observed  in  RAOB  regression  as  in  ASP  AM  regression. 
The  residual  versus  predicted  plots  for  missions  1,  4,  and  5  are  depicted  in  Figure  19. 
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(a) 


(b) 


RAOB  Pressure  (mb)  Predicted 


(c) 


RAOB  Pressure  (mb)  Predicted 


Figure  19.  Residual  versus  predicted  plots  of  RAOB  pressure  for  Global  Hawk  missions 
(a)  1,  (b)  4,  and  (c)  5. 


b.  Temperature 

All  tests  and  analyses  of  Global  Hawk  temperature  were  also  successful  and 
provided  adequate  results  to  conclude  that  Global  Hawk  temperature  reports  were  valid. 
The  linear  regression  of  Global  Hawk  temperature  with  ASPAM  and  RAOB  temperature 
produced  significant  linear  models  and  extremely  low  values  of  variance  and  RMSE. 
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1)  ASPAM 

All  five  Global  Hawk  missions  provided  772  data  points  which  were  used  in 
linear  regression  of  Global  Hawk  temperature  (x)  with  ASPAM  temperature  (y).  The 
JMP  scatterplot  and  fitted  line  is  depicted  in  Figure  20  and  the  regression  analysis  results 
are  illustrated  in  Figure  21.  The  results  include  R^,  RMSE,  ANOVA  table,  F-ratio,  and 
parameter  estimates,  along  with  the  t-test  values  for  the  hypothesis  test  =  0 .  The 

figures  show  the  Global  Hawk  is  an  accurate  predictor  of  atmospheric  temperature  as 
established  by  ASPAM. 


ASPAM  Temp  (deg  C)  Predicted  PO.OOOO 
RSq=0.99  RMSE=2.9648 


Figure  20.  Simple  linear  regression  scatterplot  and  fitted  line  of  Global  Hawk 
temperature  (x)  versus  ASPAM  temperature  (y). 


The  hypothesis  test  for  was  accomplished  using  the  JMP  output  from 
Figure  21  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
freedom  are  defined  by  n  -  2  =  770  for  Global  Hawk  versus  ASPAM  regression.  The 
JMP  values  used  were  =  0.9138359  and  =8.79  .  The  critical  t -value  and 
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were  calculated  as  n-i  -  ^0025  770  ~  ^  -963  and  =  802047.26  .  The  output  t-value 


was  to  =  -26.027  which  led  to  the  conclusion  of  rejecting  the  null  hypothesis  P^=\.  As 

previously  stated,  this  test  was  likely  to  fail  because  of  the  extremely  low  variance  of  the 
linear  model.  Therefore,  all  assumptions  of  Global  Hawk  temperature  data  quality  were 
based  on  the  regression  analysis  and  hypothesis  test  of  =  0  . 


Summary  of  Fit 


RSquare 

0.989996 

RSquare  Adj 

0.989983 

Root  Mean  Square  Error 

2.964833 

Mean  of  Response 

-31.6626 

Observations  (or  Sum  Wgts) 

772 

Analysis  of  Variance 

F  Ratio 
76196.67 
Prob  >  F 
0.0000 


Source  DF  Sum  of  Squares  Mean  Square 

Model  1  669786.45  669786 

Error  770  6768.48  9 

C.  Total  771  676554.93 


Parameter  Estimates 

Term  Estimate  Std  Error  t  Ratio  Prob>|t| 

Intercept  -5.82565  0.141941  -41.04  <.0001 

GH  Temp  (deg  C)  0.9138359  0.003311  276.04  0.0000 

Figure  21.  Regression  analysis  results  of  Global  Hawk  temperature  (x)  versus  ASPAM 
temperature  (y). 


Residual  analysis  showed  there  were  no  significant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  scale  as  the  regression  variables. 
Close  examination  of  the  data  revealed  some  possible  patterns  in  the  bias,  however. 
Global  Hawk  temperature  values  primarily  had  a  positive  bias  during  ascent  and  a 
negative  bias  during  descent.  Also,  the  amount  of  bias  seemed  to  change  in  relation  to 
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the  rate  at  whieh  the  Global  Hawk  elimbed  or  deseended.  It  is  possible  that  a  time  delay 
between  measuring  and  transmitting  atmospherie  values  by  the  Global  Hawk  would 
induee  the  bias  observed  in  this  study.  During  aseent  (deseent),  the  Global  Hawk  eould 
measure  the  temperature  at  an  altitude  lower  (higher)  than  the  height  at  whieh  it 
transmitted  the  data.  This  slight  differenee  in  flight  level  eould  lead  to  temperature 
values  higher  (lower)  than  those  at  whieh  the  data  was  relayed.  This  height  difference 
could  vary  in  relation  to  the  rate  at  which  the  Global  Hawk  climbed  or  descended  and  the 
temperature  bias  would  reflect  this  difference.  The  residual  versus  predicted  plots  for 
each  Global  Hawk  mission  are  depicted  in  Figure  22. 

2)  RAOB 

Only  three  rawinsonde  observations  were  available  for  comparison  with  Global 
Hawk  data  and  a  total  of  246  data  points  were  analyzed.  Global  Hawk  temperature  was 
used  as  a  predictor  for  RAOB  temperature  and  the  scatterplot  with  fitted  line  is  shown  in 
Figure  23.  The  JMP  regression  analysis  results  are  shown  in  Figure  24,  including  , 
RMSE,  ANOVA  table,  F-ratio,  and  parameter  estimates  with  the  t-test  results  for  the 
hypothesis  =  0  .  Regression  analysis  showed  the  Global  Hawk  is  an  accurate  predictor 
of  atmospheric  temperature  as  established  by  RAOBs. 
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ASPAM  Temp  (deg  C)  Residual  ASPAM  Temp  (deg  C)  Residual 
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ASPAM  Temp  (deg  C)  Predicted 


(d) 


ASPAM  Temp  (deg  C)  Predicted 


(e) 


ASPAM  Temp  (deg  C)  Predicted 


Figure  22.  Residual  versus  predicted  plots  of  ASPAM  temperature  for  Global  Hawk 
missions  (a-e)  1  through  5. 
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RAOB  Temp  (deg  C)  Predicted  P<.0001 
RSq=0.99  RMSE=3.1776 

Figure  23.  Simple  linear  regression  seatterplot  and  fitted  line  of  Global  Hawk 
temperature  (x)  versus  RAOB  temperature  (y). 


Summary  of  Fit 


RSquare 

0.989639 

RSquare  Adj 

0.989597 

Root  Mean  Square  Error 

3.177634 

Mean  of  Response 

-10.7955 

Observations  (or  Sum  Wgts) 

246 

Analysis  of  Variance 

Source  DF  Sum  of  Squares 

Model  1  235337.71 

Error  244  2463.76 

C.  Total  245  237801.47 


Parameter  Estimates 

Term  Estimate  Std  Error  t  Ratio  Prob>|t| 

Intercept  0.6835406  0.216101  3.16  0.0018 

GH  Temp  (deg  C)  0.9812225  0.006427  152.67  <.0001 


Mean  Square  F  Ratio 
235338  23306.86 
10  Prob  >  F 
<.0001 


Figure  24.  Regression  analysis  results  of  Global  Hawk  temperature  (x)  versus  RAOB 
temperature  (y). 


The  hypothesis  test  for  \  was  aceomplished  using  the  JMP  output  from 
Figure  24  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
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freedom  are  defined  by  n-2  =  244  for  Global  Hawk  versus  RAOB  regression.  The  JMP 
values  used  were  =0.9812225  and  =10.097.  The  eritieal  ^-value  and  SS^  were 
ealeulated  as  G/2  n-i  -  025  244  =  1  -9697  and  =  802047.26  .  The  output  t-value  was 
tg  =  -5.292  whieh  led  to  the  eonelusion  of  rejeeting  the  null  hypothesis  As 

previously  stated,  this  test  was  likely  to  fail  beeause  of  the  extremely  low  varianee  of  the 
linear  model.  Therefore,  all  assumptions  of  Global  Hawk  temperature  data  quality  were 
based  on  the  regression  analysis  and  hypothesis  test  of  =  0  . 

Residual  analysis  showed  there  were  no  signifieant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  seale  as  the  regression  variables. 
Similar  patterns  in  the  bias  were  observed  in  RAOB  regression  as  in  ASP  AM  regression. 
The  residual  versus  predieted  plots  for  missions  1,  4,  and  5  are  depieted  in  Figure  25. 

c.  Wind  speed 

All  tests  and  analyses  of  Global  Hawk  wind  speed  were  also  sueeessful  and 
provided  adequate  results  to  eonelude  that  Global  Hawk  wind  speed  reports  were  valid. 
The  linear  regression  of  Global  Hawk  wind  speed  with  ASPAM  and  RAOB  wind  speed 
produced  significant  linear  models  and  extremely  low  values  of  variance  and  RMSE. 
Several  studies  have  compared  AGARS  wind  data  with  other  sources,  such  as 
rawinsondes,  and  the  amount  of  error  ranged  from  1.1  m  s’'  to  3.1  m  s‘'  (Schwartz  and 
Benjamin  1995;  Richner  and  Gutermann  1987;  Benjamin  et  al.  1999;  Morone  1986;  Nash 
1994).  Nash  (1994)  attributed  1-1.5  m  s'^  of  the  total  2.6  m  s"'  wind  error  to  instrument 
uncertainty  and  the  rest  to  atmospheric  variability  between  the  locations  of  the  compared 
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measurements.  The  Global  Hawk  wind  error  in  the  following  tests  is  within  this  range  of 
variability  and,  therefore,  is  assumed  to  be  as  accurate  and  as  useful  as  ACARS  data 
currently  in  use. 
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Figure  25.  Residual  versus  predicted  plots  of  RAOB  temperature  for  Global  Hawk 
missions  (a)  1,  (b)  4,  and  (c)  5. 


1)  ASPAM 

All  five  Global  Hawk  missions  provided  772  data  points  which  were  used  in 
linear  regression  of  Global  Hawk  wind  speed  (x)  with  ASPAM  wind  speed  (y).  The  JMP 
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scatterplot  and  fitted  line  is  depicted  in  Figure  26  and  the  regression  analysis  results  are 
illustrated  in  Figure  27.  The  results  include  R^,  RMSE,  ANOVA  table,  F-ratio,  and 
parameter  estimates,  along  with  the  t-test  values  for  the  hypothesis  test  =  0  .  The 


figures  show  the  Global  Hawk  is  an  accurate  predictor  of  atmospheric  wind  speed  as 
established  by  ASP  AM. 


ASPAM  Wind  Speed  (m/s)  Predicted 
P<.0001  RSq=0.82  RMSE=2.9589 

Figure  26.  Simple  linear  regression  scatterplot  and  fitted  line  of  Global  Hawk  wind 
speed  (x)  versus  ASPAM  wind  speed  (y). 

The  hypothesis  test  for  /5^=  \  was  accomplished  using  the  JMP  output  from 
Figure  27  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
freedom  are  defined  by  n  -  2  =  770  for  Global  Hawk  versus  ASPAM  regression.  The 
JMP  values  used  were  =  0.5581524  and  =  33.053  .  The  critical  /-value  and 
were  calculated  as  «-2  =  hais  770  =  1  -963  and  =  95467.40  .  The  output  /-value  was 
/g  =  -23.746  which  led  to  the  conclusion  of  rejecting  the  null  hypothesis  =  1 .  As 
previously  stated,  this  test  was  likely  to  fail  because  of  the  extremely  low  variance  of  the 
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linear  model.  Therefore,  all  assumptions  of  Global  Hawk  wind  speed  data  quality  were 
based  on  the  regression  analysis  and  hypothesis  test  of  =  0  . 


■ 

Summaty  of  Fit 


RSquare 

0.815215 

RSquare  Adj 

0.814975 

Root  Mean  Square  Error 

2.958914 

Mean  of  Response 

8.96237 

Observations  (or  Sum  Wgts) 

772 

Analysis  of  Variance 

F  Ratio 
3397.004 
Prob  >  F 
<•0001 


Source  DF  Sum  of  Squares  Mean  Square 

Model  1  29741.353  29741.4 

Error  770  6741.483  8.8 

C.  Total  771  36482.836 


■  ■ 

Parameter  Estimates 

Term  Estimate  Std  Error  t  Ratio  Prob>|t| 

Intercept  0.3511546  0.182126  1.93  0.0542 

GH  Wind  Speed  (m/s)  0.5581524  0.009576  58.28  <.0001 

Figure  27.  Regression  analysis  results  of  Global  Hawk  wind  speed  (x)  versus  ASPAM 
wind  speed  (y). 


Residual  analysis  showed  there  were  no  signifieant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  scale  as  the  regression  variables. 
Close  examination  of  the  data  revealed  some  possible  patterns  in  the  bias,  however.  A 
majority  of  the  extreme  bias  values  ( -2  m  5  '  >  bias  >  2  m  5  ' )  were  reported  in  the  lower 
( <  1 0  Kft )  and  upper  ( >  40  Kft )  levels  of  the  Global  Hawk  missions.  It  is  the  assumption 
of  this  research  that  the  bias  was  induced  by  the  mesoscale  variability  of  the  winds  in  the 
lower  levels  and  significantly  higher  wind  speeds  in  the  upper  levels.  Global  Hawk  wind 
measurements  could  be  affected  by  higher  upper-level  wind  speeds  and  these  impacts 
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would  vary  with  the  speed  and  attitude  of  the  aireraft  relative  to  the  mean  flow.  An 
additional  eonsideration  is  rawinsonde  balloon  drift  in  the  upper-levels,  whieh  would 
influenoe  ASPAM  profdes.  The  residual  versus  predieted  plots  are  depicted  in  Figure  28. 
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Figure  28.  Residual  versus  predicted  plots  of  ASPAM  wind  speed  for  Global  Hawk 
missions  (a-e)  1  through  5 . 
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2)  RAOB 


Only  three  rawinsonde  observations  were  available  for  eomparison  with  Global 
Hawk  data  and  a  total  of  246  data  points  were  analyzed.  Global  Hawk  wind  speed  was 
used  as  a  predietor  for  RAOB  wind  speed  and  the  seatterplot  with  fitted  line  is  shown  in 
Figure  29.  The  JMP  regression  analysis  results  are  shown  in  Figure  30,  ineluding  R^, 
RMSE,  ANOVA  table,  F-ratio,  and  parameter  estimates  ineluding  the  t-test  results  for  the 
hypothesis  =  0  .  Regression  analysis  results  show  the  Global  Hawk  is  an  accurate 
predictor  of  atmospheric  wind  speed  as  established  by  RAOBs. 


RAOB  Wind  Speed  (m/s)  Predicted 
P<.0001  RSq=0.80  RMSE=3.2097 

Figure  29.  Simple  linear  regression  seatterplot  and  fitted  line  of  Global  Hawk  wind 
speed  (x)  versus  RAOB  wind  speed  (y). 


The  hypothesis  test  for  \  was  accomplished  using  the  JMP  output  from 
Figure  24  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
freedom  are  defined  by  n-2  =  244  for  Global  Hawk  versus  RAOB  regression.  The  JMP 
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values  used  were  =0.5366051  and  a  =10.30.  The  eritieal^-value  and  SS^  were 
calculated  as  «-2  “  ^0025  770  ~  1  -9697  and  SS^  =  95467.40  .  The  output  t-value  was 
tg  =  -44.613  which  led  to  the  conclusion  of  rejecting  the  null  hypothesis  =  1 .  As 

previously  stated,  this  test  was  likely  to  fail  because  of  the  extremely  low  variance  of  the 
linear  model.  Therefore,  all  assumptions  of  Global  Hawk  temperature  data  quality  were 
based  on  the  regression  analysis  and  hypothesis  test  of  =  0  . 


Summary  of  Fit 


RSquare 

0.796424 

RSquare  Adj 

0.795589 

Root  Mean  Square  Error 

3.209707 

Mean  of  Response 

9.1 

Observations  (or  Sum  Wgts) 

246 

Analysis  of  Variance 

F  Ratio 
954.5672 
Prob  >  F 
<.0001 


Source  DF  Sum  of  Squares  Mean  Square 

Model  1  9834.163  9834.16 

Error  244  2513.742  10.30 

C.  Total  245  12347.904 


Parameter  Estimates 


Term 

Intercept 

GH  Wind  Speed  (m/s) 


Estimate  Std  Error 
1.6961335  0.315127 
0.5366051  0.017368 


t  Ratio  Prob>|t| 
5.38  <.0001 

30.90  <.0001 


Figure  30.  Regression  analysis  results  of  Global  Hawk  wind  speed  (x)  versus  RAOB 
wind  speed  (y). 


Residual  analysis  showed  there  were  no  significant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  scale  as  the  regression  variables. 
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Similar  patterns  in  the  bias  were  observed  in  RAOB  regression  as  in  ASP  AM  regression. 
The  residual  versus  predicted  plots  for  missions  1,  4,  and  5  are  depicted  in  Figure  31. 
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Figure  31.  Residual  versus  predicted  plots  of  RAOB  wind  speed  for  Global  Flawk 
missions  (a)  1,  (b)  4,  and  (c)  5. 


3)  ASPAM  -  ADJUSTED  WIND  DATA 

Investigation  of  the  data  revealed  certain  wind  data  points  which  could  be 
excluded  from  analysis  based  on  certain  assumptions.  The  extreme  outliers  in  the 
original  regression  were  excluded  due  to  the  observation  that  they  were  grouped  together 
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during  specific  periods  of  each  of  the  missions.  It  was  assumed  that  either  the  data  were 
relayed  incorrectly  due  to  a  system  problem  or  that  an  unknown  atmospheric  occurrence 
corrupted  the  information.  A  total  of  122  extreme  outliers  were  excluded  based  on  these 
assumptions. 

Wind  measurements  could  be  degraded  due  to  the  motion  and  flight  attitude  of  the 
measuring  airframe  (Axford  1968).  Painting  (2002)  determined  an  uncertainty  of  wind 
data  to  be  2-3  m  s'\  which  increased  during  aircraft  maneuvers.  Bisiaux  et  al.  (1983) 
also  calculated  higher  errors  in  wind  data  during  aircraft  maneuvers  and  excluded  reports 
with  roll  angles  greater  than  5  degrees.  This  research  excluded  data  points  at  which  the 
aircraft  was  either  pitching  or  rolling  in  excess  of  10  degrees  based  on  the  assumptions  of 
aircraft  maneuvering  impacts  on  wind  data.  A  total  of  129  reports  were  excluded  for 
extreme  pitch  and  254  reports  were  excluded  for  extreme  roll  of  the  Global  Hawk. 

Data  points  were  also  excluded  based  on  the  assumptions  of  erroneous  Global 
Hawk  data  as  described  in  Section  4.b.l.  The  system  reported  to  be  in  Altitude  Hold 
Mode  but  the  aircraft  altitude  continued  to  fluctuate,  possibly  indicating  suspect  data. 
Notes  pertaining  to  the  Global  Hawk  data  and  varying  flight  modes  are  shown  in 
Appendix  D.  A  total  of  68  reports  were  excluded  based  on  these  assumptions. 

In  all,  440  data  points  were  excluded  from  statistical  analysis  based  on  these 
assumptions.  Some  of  the  reports  fell  into  two  or  more  of  the  categories  listed  above, 
leading  to  a  lower  number  of  exclusions  than  expected.  Of  particular  interest  was  that  all 
but  one  data  point  from  mission  1  were  excluded.  The  regression  analysis  and  hypothesis 
test  results  follow. 
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All  five  Global  Hawk  missions  provided  332  data  points  wtiieli  were  used  in 
linear  regression  of  Global  Hawk  wind  speed  (x)  with  ASPAM  wind  speed  (y).  The  JMP 
seatterplot  and  fitted  line  is  depieted  in  Figure  32  and  the  regression  analysis  results  are 
illustrated  in  Figure  33.  The  results  inelude  R^,  RMSE,  ANOVA  table,  F-ratio,  and 
parameter  estimates,  along  with  the  ^-test  values  for  the  hypothesis  test  =  0  .  The 

figures  show  the  Global  Hawk  is  an  aeeurate  predietor  of  atmospherie  wind  speed  as 
established  by  ASPAM. 


ASPAM  Wind  Speed  (m/s)  Predicted 
P<.0001  RSq=0.90  RMSE=2.3006 


Figure  32.  Simple  linear  regression  seatterplot  and  fitted  line  of  Global  Hawk  wind 
speed  (x)  versus  ASPAM  wind  speed  (y)  with  adjusted  wind  data. 


The  hypothesis  test  for  was  accomplished  using  the  JMP  output  from 

Figure  33  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
freedom  are  defined  by  n  -  2  =  330  for  Global  Hawk  versus  ASPAM  regression.  The 
JMP  values  used  were  =  0.5969934  and  =  5.29 .  The  critical  /-value  and  SS^ 
were  calculated  as  «-2  =  025  330  ~  ^  -9672  and  SS^  =  45653 .76  .  The  output  /-value 
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was  =  -37.439  which  led  to  the  eonelusion  of  rejeeting  the  null  hypothesis  =  1 .  As 

previously  stated,  this  test  was  likely  to  fail  because  of  the  extremely  low  variance  of  the 
linear  model.  Therefore,  all  assumptions  of  Global  Hawk  wind  speed  data  quality  were 
based  on  the  regression  analysis  and  hypothesis  test  of  =  0  . 


Summary  of  Fit 


RSquare 

0.903061 

RSquare  Adj 

0.902768 

Root  Mean  Square  Error 

2.300597 

Mean  of  Response 

10.35633 

Observations  (or  Sum  Wgts) 

332 

Analysis  of  Variance 

Source 

DF  Sum  of  Squares 

Mean  Square  F  Ratio 

Model 

1 

16271.052 

16271.1  3074.217 

Error 

330 

1746.607 

5.3  Prob  >  F 

C.  Total 

331 

18017.658 

<•0001 

Parameter  Estimates 

Term 

Estimate 

Std  Error 

t  Ratio  Prob>|t 

Intercept 

0.5407221 

0.217444 

2.49  0.0134 

GH  Wind  Speed  (m/s) 

0.5969934 

0.010767 

55.45  <.0001 

Figure  33.  Regression  analysis  results  of  Global  Hawk  wind  speed  (x)  versus  ASPAM 
wind  speed  (y)  with  adjusted  wind  data. 


Residual  analysis  showed  there  were  no  significant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  scale  as  the  regression  variables. 
Close  examination  of  the  data  revealed  some  possible  patterns  in  the  bias,  however.  A 
majority  of  the  extreme  bias  values  ( -2  m  5  '  >  bias  >  2  m  5  ' )  were  reported  in  the  lower 
( <  1 0  Kft )  and  upper  (  >  40  Kft )  levels  of  the  Global  Hawk  missions.  It  is  the  assumption 
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of  this  research  that  the  bias  was  induced  by  the  mesoscale  variability  of  the  winds  in  the 
lower  levels  and  significantly  higher  wind  speeds  in  the  upper  levels.  Global  Hawk  wind 
measurements  could  be  affected  by  higher  upper-level  wind  speeds  and  these  impacts 
would  vary  with  the  speed  and  attitude  of  the  aircraft  relative  to  the  mean  flow.  An 
additional  consideration  is  rawinsonde  balloon  drift  in  the  upper-levels,  which  would 
influence  ASPAM  profiles.  The  residual  versus  predicted  plots  for  each  Global  Hawk 
mission  are  depicted  in  Figure  34.  Only  one  data  point  was  available  for  mission  1  after 
excluding  suspect  wind  data.  Therefore,  the  residual  plot  for  mission  1  is  empty. 

4)  RAOB  -  ADJUSTED  WIND  DATA 

The  same  assumptions  were  made  and  rows  excluded  as  in  Section  5.C.3.  A  total 
of  172  data  rows  were  excluded  leaving  only  74  rows  of  data  for  analysis  with 
rawinsonde  observations.  Global  Hawk  wind  speed  was  used  as  a  predictor  for  RAOB 
wind  speed  and  the  scatterplot  with  fitted  line  is  shown  in  Figure  35.  The  JMP  regression 
analysis  results  are  shown  in  Figure  36,  including  ,  RMSE,  ANOVA  table,  F-ratio,  and 
parameter  estimates  including  the  t-test  results  for  the  hypothesis  =  0  .  Regression 

analysis  results  show  the  Global  Hawk  is  an  accurate  predictor  of  atmospheric  wind 
speed  as  established  by  RAOBs. 
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Figure  34.  Residual  versus  predicted  plots  of  ASPAM  wind  speed  for  Global  Hawk 
missions  (a-e)  1  through  5  with  adjusted  wind  data. 


64 


The  hypothesis  test  for  \  was  accomplished  using  the  JMP  output  from 

Figure  36  and  calculations  in  Mathcad  software.  Recall  that  a  =  0.05  and  the  degrees  of 
freedom  are  defined  by  n-2  =  12  for  Global  Hawk  versus  RAOB  regression  with 
adjusted  wind  data.  The  JMP  values  used  were  =  0.6029932  and  <j^  =  4.37 .  The 
critical  t-value  and  SS^  were  calculated  as  =^002572  =1-9935  and 

=  45653.76 .  The  output  t-value  was  tg  =  -40.578  which  led  to  the  conclusion  of 
rejecting  the  null  hypothesis  =  1 .  As  previously  stated,  this  test  was  likely  to  fail 

because  of  the  extremely  low  variance  of  the  linear  model.  Therefore,  all  assumptions  of 
Global  Hawk  temperature  data  quality  were  based  on  the  regression  analysis  and 
hypothesis  test  of  =  0  . 


0  5  10  15  20  25 

RAOB  Wind  Speed  (m/s)  Predicted 
P<.0001  RSq=0.95  RMSE=2.0899 


Figure  35.  Simple  linear  regression  scatterplot  and  fitted  line  of  Global  Hawk  wind 
speed  (x)  versus  RAOB  wind  speed  (y)  with  adjusted  wind  data. 
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Summary  of  Fit 


RSquare 
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RSquare  Adj 
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Root  Mean  Square  Error 
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Mean  of  Response 
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Observations  (or  Sum  Wgts) 

74 

Analysis  of  Variance 

F  Ratio 
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Prob  >  F 
<.0001 


Source  DF  Sum  of  Squares  Mean  Square 

Model  1  5989.5106  5989.51 

Error  72  314.4798  4.37 

C.  Total  73  6303.9905 


Parameter  Estimates 


Term 

Intercept 

GH  Wind  Speed  (m/s) 


Estimate  Std  Error 
0.1362796  0.361085 
0.6029932  0.016283 


t  Ratio  Prob>|t| 
0.38  0.7070 

37.03  <.0001 


Figure  36.  Regression  analysis  results  of  Global  Flawk  wind  speed  (x)  versus  RAOB 
wind  speed  (y)  with  adjusted  wind  data. 


Residual  analysis  showed  there  were  no  significant  problems  in  the  regression. 
The  residuals  were  small  when  examined  on  the  same  scale  as  the  regression  variables. 
Similar  patterns  in  the  bias  were  observed  in  RAOB  regression  as  in  ASP  AM  regression. 
The  residual  versus  predicted  plots  for  missions  1,  4,  and  5  are  depicted  in  Figure  37. 
Only  one  data  point  was  available  for  mission  1  after  excluding  suspect  wind  data. 
Therefore,  the  residual  plot  for  mission  1  is  empty. 
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Figure  37.  Residual  versus  predicted  plots  of  RAOB  wind  speed  for  Global  Flawk 
missions  (a)  1,  (b)  4,  and  (c)  5  with  adjusted  wind  data. 


d.  Wind  direction 

1)  RMS  VECTOR  ERROR 

The  root  mean  square  vector  error  of  the  Global  Flawk  winds  versus  the  ASPAM 
and  RAOB  winds  was  calculated  using  both  the  original  and  the  adjusted  datasets.  The 
wind  direction  and  wind  speed  were  taken  into  account  using  this  method.  The  RMS 
vector  error  of  Global  Hawk  versus  ASPAM  winds  was  0.806  m  s'*  with  the  original 
dataset  and  0.566  m  s'*  using  the  adjusted  data.  The  RMS  vector  error  of  Global  Hawk 
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versus  RAOB  winds  was  1.403  m  s'*  with  the  original  dataset  and  0.988  m  s'*  using  the 
adjusted  data. 


2)  Mean  DIFFERENCE 

The  ASPAM  and  RAOB  wind  direetion  values  were  subtraeted  from  the  Global 
Hawk  wind  direetion  measurements  using  the  teehnique  illustrated  in  Seetion  4.e.2.  The 
mean  differenees  between  Global  Hawk  and  ASPAM  data  were  3 1 .4°  for  the  original 
dataset  and  26.3°  for  the  adjusted  data.  The  mean  differenees  between  the  Global  Hawk 
and  RAOB  data  were  36.5°  for  the  original  dataset  and  42.4°  for  the  adjusted  data.  The 
degradation  in  the  RAOB  differenee  with  the  adjusted  dataset  was  most  likely  due  to  the 
signifieantly  lower  number  of  observations  available  after  the  suspeet  information  was 
eliminated. 

A  signifieant  portion  of  the  wind  direetion  error  ean  be  attributed  to  light  wind 
eonditions,  whieh  were  determined  as  5  m  s'*  or  less  for  this  researeh.  Light  winds 
oeeurred  about  64%  of  the  time  with  ASPAM  and  RAOB  reports  and  14%  of  the  time  for 
Global  Hawk  reports.  Sehwartz  and  Benjamin  (1995)  determined  a  mean  direetional 
differenee  of  35.1°  when  eomparing  AGARS  and  rawinsonde  wind  data,  whieh  is 
eomparable  to  the  results  of  this  researeh. 
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6.  Conclusions  and  Recommendations 


a.  Conclusions 

The  Global  Hawk  Air  Dara  System  is  an  aeeurate  predietor  of  the  atmospherie 
variables  of  pressure,  temperature,  wind  speed,  and  wind  direetion.  The  statistieal  tests 
performed  in  this  researeh  show  an  almost  perfeet  fit  of  Global  Hawk  pressure  and 
temperature  data  with  both  ASPAM  and  RAOB  values.  The  linear  fit  for  wind  data  was 
also  signifieant,  espeeially  after  some  assumptions  and  exelusions  were  made.  The 
results  demonstrate  environmental  intelligenee  gathered  by  the  Global  Hawk  is  preeise 
and,  therefore,  of  value  to  both  the  intelligenee  and  meteorologieal  eommunities. 

The  linear  fit  of  Global  Hawk  pressure  with  ASPAM  and  RAOB  data  was 
virtually  exaet  with  values  of  0.999.  The  RMSE  values  were  exeeptionally  low  with 
8.29  mb  for  ASPAM  regression  and  5.52  mb  for  RAOB  regression. 

Simple  linear  regression  revealed  Global  Hawk  temperature  eomparisons  with 
ASPAM  and  RAOB  reports  were  also  very  preeise  with  values  of  0.990.  The  RMSE 
values  were  also  small  with  2.96  °C  for  ASPAM  regression  and  3.18  °C  for  RAOB 
regression. 

The  statistieal  analysis  of  Global  Hawk  wind  speed  data  was  informative  and  the 
outeome  was  eomparable  to  aeeuraey  tests  performed  on  AGARS  data.  Initial  regression 
with  ASPAM  and  RAOB  data  provided  R^  values  of  0.815  and  0.796,  respeetively.  The 
RMSE  values  were  2.96  m  s"'  for  ASPAM  regression  and  3.21  m  s'^  for  RAOB 
regression.  After  adjusting  the  data  and  making  exelusions  based  on  eertain  assumptions, 
the  results  improved  drastieally.  The  exeluded  data  eonsisted  of  extreme  outliers  that 
were  grouped  together,  exeessive  piteh  and  roll  of  the  Global  Hawk,  and  possible 
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erroneous  Global  Hawk  reports.  The  new  values  were  0.903  for  ASP  AM  and  0.950 
for  RAOB.  The  adjusted  RMSE  values  for  ASPAM  and  RAOB  regression  were  2.30  m 
s'^  and  2.09  m  s'\  respeetively.  This  typieal  error  magnitude  of  the  Global  Hawk  wind 
speed  data  is  similar  to  the  results  of  numerous  quality  assessments  of  AGARS  data. 

Wind  direetion  data  eould  not  be  analyzed  by  normal  methods  due  to  its  non¬ 
linearity  but  useful  information  was  obtained  by  ealeulating  the  RMS  veetor  errors  and 
mean  differenees.  Analysis  with  the  original  datasets  produeed  RMS  veetor  errors  of 
0.806  m  s'^  for  Global  Hawk  versus  ASPAM  and  1.403  m  s"'  for  Global  Hawk  versus 
RAOB.  The  mean  differenees  in  direetion  between  the  initial  Global  Hawk  and  “ground 
truth”  datasets  were  31.4°  with  ASPAM  and  36.5°  with  RAOB.  After  assumptions  were 
made  and  the  data  were  adjusted,  RMS  veetor  errors  improved  to  0.566  m  s'^  for  Global 
Hawk  versus  ASPAM  and  0.988  m  s'^  for  Global  Hawk  versus  RAOB.  The  adjusted  data 
produeed  different  results  for  the  mean  differenee  in  wind  direetion.  The  mean  differenee 
improved  to  26.3°  for  the  ASPAM  data  but  worsened  to  42.4°  for  the  RAOB  data. 
Independent  verifieation  of  AGARS  wind  data  has  produeed  RMS  veetor  errors  ranging 
from  1.1  m  s'*  to  3.1  m  s'*  and  mean  differenees  between  30°  and  40°.  Mean  differenee 
values  this  large  should  raise  some  doubt  as  to  the  aeeuraey  of  the  wind  data.  An 
assumption  was  made  that  light  winds  (<  5  m  s'*)  would  introduee  variability  in  the  wind 
direetion  data.  RMS  veetor  error  ealeulations  alleviate  any  problems  due  to  light  and 
variable  winds  beeause  both  wind  direetion  and  wind  speed  are  taken  into  aeeount.  The 
Global  Hawk  RMS  veetor  errors  produeed  in  this  researeh  are  analogous  with  those  of 
AGARS  data  found  in  previous  studies  (Sehwartz  and  Benjamin  1995;  Riehner  and 
Gutermann  1987;  Benjamin  et  al.  1999;  Morone  1986;  Nash  1994).  Even  so,  the 
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accuracy  of  ACARS  wind  information  has  been  deemed  adequate  to  suoeessfully 
enhanee  foreeast  produets  and  NWP  models  for  the  past  deeade. 

The  results  above  demonstrate  the  aeeuraey  of  the  Global  Hawk  environmental 
information.  Several  assumptions  were  made  in  this  researeh  but  the  overall  results 
provided  suffieient  justifieation  that  this  UAV  would  serve  well  as  a  meteorologieal 
sensing  platform.  Alterations  to  the  Global  Hawk  aireraft  eonfiguration  or  mission 
profile  are  not  required  in  order  to  obtain  and  exploit  this  data.  Therefore,  steps  should 
be  taken  to  ensure  this  souree  of  environmental  information  is  made  available  to  military 
weather  ageneies  as  well  as  UAV  operators  and  planners.  An  inerease  in  aeeurate 
weather  data  would  lead  to  improved  mission  effeetiveness,  espeeially  over  data-sparse 
regions.  The  Global  Hawk  has  the  ability  to  provide  atmospherie  soundings  during  its 
aseent  and  deseent  and  its  near  real  time  data  eould  benefit  mission  operators,  planners, 
and  weather  foreeasters  alike.  Sample  soundings  from  all  five  Global  Hawk  missions  are 
depieted  in  Figure  38,  illustrating  just  one  of  the  many  possible  benefits  of  aequiring  and 
utilizing  this  aeeurate  and  timely  information.  Dual  lines  on  some  of  the  plots  indieate 
two  soundings  were  ereated,  one  during  aseent  and  the  other  during  deseent.  The 
tropopause,  stratospherie  inversion,  and  low-level  inversions  are  evident  on  several  plots. 
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b.  Recommendations 


The  RQ-4A  Global  Hawk  is  a  proven  asset  to  the  intelligenee  and  operational 
eommunities  and  eould  benefit  the  meteorologieal  eommunity  as  well,  if  the  atmospherie 
data  from  the  system  were  made  available.  Weather  data  are  measured  and  transmitted 
regularly  by  the  Global  Hawk  during  all  phases  of  its  operation  so  efforts  must  be  made 
to  aequire  this  readily  abundant  information.  Additional  sensors  or  meteorologieal 
payloads  are  not  required  to  make  use  of  this  valuable  data.  Thus,  it  appears  the  only 
thing  preventing  the  use  of  this  atmospherie  intelligenee  is  eoordination  between  the 
affeeted  ageneies  to  determine  an  aequisition  and  eommunieations  systems  strategy, 
while  maintaining  an  awareness  of  the  elassified  nature  of  the  majority  of  Global  Hawk 
flight  routes  and  missions. 

1)  For  Air  Force  Weather 

Steps  should  be  taken  to  obtain  Global  Hawk  environmental  intelligenee  in  near 
real  time  from  the  mission  eontrol  elements  or  the  airframe  itself  Coordination  with  the 
Global  Hawk  System  Program  Offiee  and  the  88*  Weather  Squadron  at  Wright-Patterson 
AFB  is  paramount  to  ensure  aequisition  of  this  valuable  weather  intelligenee,  espeeially 
in  data-sparse  battlefield  environments.  Atmospherie  information  are  available  in 
numerous  data-sparse  regions  during  Global  Hawk  operations  and  this  intelligenee  eould 
be  used  to  enhanee  model  initialization,  improve  foreeast  produets,  and  provide 
additional  data  points  for  analysis. 
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2)  For  Global  Hawk  SPO 

Correspond  with  Air  Force  weather  agencies  to  make  Global  Hawk  weather  data 
available  to  meteorologists  and  mission  planners.  Doing  so  would  not  only  improve  the 
weather  support  provided  to  Global  Hawk  operators  but  also  allow  trained  technicians  to 
evaluate  the  atmospheric  reports  and  ensure  operations  are  completed  successfully. 

Analyze  the  datasets  used  in  this  research  along  with  data  from  other  missions  to 
possibly  determine  the  causes  of  the  bias  in  the  data.  Research  the  information  provided 
in  Appendix  D  to  resolve  any  possible  measurement  or  transmission  problems  within  the 
Global  Hawk  system. 

3)  For  FUTURE  RESEARCH 

Several  additional  tests  and  techniques  were  contemplated  during  this  research  but 
time  was  a  major  impact  not  only  on  the  methodology  but  also  the  final  thesis  product. 
More  intense  regression  analysis  of  all  datasets  should  be  explored  as  well  as  other 
possible  statistical  tests.  Further  investigation  into  circular  data  analysis  could  provide 
useful  information  as  well  as  a  new  approach  to  verifying  environmental  data. 

Colleagues  and  professors  had  very  interesting  ideas  and  the  question,  “Have  you 
considered. . .?”  was  asked  on  several  occasions.  This  meteorologist  leaves  many 
considerations  and  unanswered  questions  to  future  students  and  fellow  researchers. 
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Appendix  A;  Acronyms 


ACARS 

Aircraft  Communications  Addressing  and  Reporting  System 

ADS 

Air  Data  System 

AFB 

Air  Eorce  Base 

AFCCC 

Air  Eorce  Combat  Climatology  Center 

AFIT 

Air  Eorce  Institute  of  Technology 

AIREPS 

Alcraft  REPortS 

AMDAR 

Aircraft  Meteorological  DAta  Relay 

ANOVA 

Analysis  of  Variance 

ARINC 

Aeronautical  Radio,  Inc 

ASPAM 

Atmospheric  Slant  Path  Analysis  Model 

AV7 

Air  Vehicle  7 

COAMPS 

Coupled  Ocean/ Atmospheric  Mesoscale  Prediction  System 

CONUS 

CONtinental  United  States 

CWSU 

Center  Weather  Service  Unit 

DEVG 

Derived  Equivalent  Vertical  Gust 

DoD 

Department  of  Defense 

DEED 

Digital  Terrain  Elevation  Data 

EDR 

Eddy  Dissipation  Rate 

EETS 

EEight  Test  Squadron 

ENMOC 

Elect  Numerical  Meteorology  and  Oceanography  Center 

ESE 

Eorecast  Systems  Eaboratory 

IMMC 

Integrated  Mission  Management  Computer 

MSE 

Mean  Square  Error 

MSE 

Mean  Sea  Eevel 

MVOI 

Multivariate  Optimal  Interpolation 

NOAA 

National  Oceanic  and  Atmospheric  Administration 

NOGAPS 

Navy  Operational  Global  Atmospheric  Prediction  System 

NWP 

Numerical  Weather  Prediction 

NWSEO 

National  Weather  Service  Eorecast  Offices 

PIREPS 

Pilot  REPortS 

RAOB 

RAwinsonde  OBservation 

RMSE 

Root  Mean  Square  Error 

RTNEPH 

Real-Time  NEPHanalysis 

SPC 

Storm  Prediction  Center 

SPO 

System  Program  Office 

TDY 

Temporary  DutY 

WPAEB 

Wright-Patterson  Air  Eorce  Base 

WRE 

Weather  Research  and  Eorecasting  model 

UAV 

Unmanned  Aerial  Vehicle 

US 

United  States 

UTC 

Universal  Time  Code 
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Appendix  B:  Upper-Level  Wind  Analyses 


Following  are  upper-level  wind  analyses  corresponding  with  the  five  Global 
Hawk  missions  used  in  this  research.  Charts  are  centered  over  California,  represent  both 
300-mb  and  200-mb  analyses,  and  are  valid  at  1200  and  0000  UTC  on  the  dates  of  the 
Global  Hawk  flights  (National  Weather  Service  Storm  Prediction  Center  2003). 


Figure  B-1.  Upper-level  wind  analyses  for  29  -  30  May  2003.  (a)  300  mb  at  29/1200 
UTC,  (b)  300  mb  at  30/0000  UTC,  (c)  250  mb  at  29/1200  UTC,  (d)  250  mb  at  30/0000 
UTC  (National  Weather  Service  Storm  Prediction  Center  2003). 
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(a)  (b) 


Figure  B-2.  Upper-level  wind  analyses  for  03  -  04  June  2003.  (a)  300  mb  at  03/1200 
UTC,  (b)  300  mb  at  04/0000  UTC,  (c)  250  mb  at  03/1200  UTC,  (d)  250  mb  at  04/0000 
UTC  (National  Weather  Service  Storm  Prediction  Center  2003). 
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Figure  B-3.  Upper-level  wind  analyses  for  26  -  27  June  2003.  (a)  300  mb  at  26/1200 
UTC,  (b)  300  mb  at  27/0000  UTC,  (c)  250  mb  at  26/1200  UTC,  (d)  250  mb  at  27/0000 
UTC  (National  Weather  Service  Storm  Prediction  Center  2003). 
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(a)  (b) 


Figure  B-4.  Upper-level  wind  analyses  for  08  -  09  August  2003.  (a)  300  mb  at 
08/1200  UTC,  (b)  300  mb  at  09/0000  UTC,  (c)  250  mb  at  08/1200  UTC,  (d)  250  mb  at 
09/0000  UTC  (National  Weather  Service  Storm  Prediction  Center  2003). 
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Figure  B-5.  Upper-level  wind  analyses  for  15  -  16  August  2003.  (a)  300  mb  at 
15/1200  UTC,  (b)  300  mb  at  16/0000  UTC,  (c)  250  mb  at  15/1200  UTC,  (d)  250  mb  at 
16/0000  UTC  (National  Weather  Service  Storm  Prediction  Center  2003). 
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Appendix  C:  Global  Hawk  Flight  Level  Information 


Following  are  the  Global  Hawk  flight  levels  used  in  this  study  and  the  differenee 
between  the  aetual  Global  Hawk  altitudes  and  the  eorresponding  500-foot  levels. 
Highlighted  rows  denote  the  flight  levels  with  the  maximum  differenee  for  the  aseending 
and  deseending  segments  of  eaeh  mission. 


Global  Hawk 
Mission  # 

Global  Hawk 
Flight  Mode 

Global  Hawk 
Flight  Level 
(ft  MSL) 

500-Foot 
Flight  Level 
(ft  MSL) 

Difference 
(ft  MSL) 

1 

Aseent 

2507.35 

2500 

7.35 

1 

Aseent 

2985.78 

3000 

14.22 

1 

Aseent 

3507.04 

3500 

7.04 

1 

Aseent 

3984.49 

4000 

15.51 

1 

Aseent 

4514.88 

4500 

14.88 

1 

Aseent 

5002.88 

5000 

2.88 

1 

Aseent 

5518.91 

5500 

18.91 

1 

Aseent 

5988.93 

6000 

11.07 

1 

Aseent 

6504.54 

6500 

4.54 

1 

Aseent 

7000.45 

7000 

0.45 

1 

Aseent 

7512.73 

7500 

12.73 

1 

Aseent 

8024.41 

8000 

24.41 

1 

Aseent 

8499.48 

8500 

0.52 

1 

Aseent 

8981.06 

9000 

18.94 

1 

Aseent 

9518.07 

9500 

18.07 

1 

Aseent 

10010.34 

10000 

10.34 

1 

Aseent 

10518.62 

10500 

18.62 

1 

Aseent 

10981.53 

11000 

18.47 

1 

Aseent 

11503.16 

11500 

3.16 

1 

Aseent 

12023.56 

12000 

23.56 

1 

Aseent 

12480.17 

12500 

19.83 

1 

Aseent 

12974.53 

13000 

25.47 

1 

Aseent 

13516.05 

13500 

16.05 

1 

Aseent 

14006.66 

14000 

6.66 

1 

Aseent 

14479.17 

14500 

20.83 

1 

Aseent 

15017.81 

15000 

17.81 

1 

Aseent 

15501.77 

15500 

1.77 

81 


I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 


Ascent 

15972.29 

16000 

27.71 

Ascent 

16498.79 

16500 

1.21 

Aseent 

16983.57 

17000 

16.43 

Aseent 

17474.67 

17500 

25.33 

Ascent 

17990.14 

18000 

9.86 

Ascent 

18503.91 

18500 

3.91 

Ascent 

19002.34 

19000 

2.34 

Aseent 

19517.82 

19500 

17.82 

Ascent 

19993.80 

20000 

6.20 

Descent 

19988.59 

20000 

11.41 

Descent 

19505.87 

19500 

5.87 

Descent 

18991.44 

19000 

8.56 

Descent 

18528.68 

18500 

28.68 

Descent 

17996.63 

18000 

3.37 

Descent 

17512.96 

17500 

12.96 

Descent 

16986.71 

17000 

13.29 

Descent 

16502.65 

16500 

2.65 

Descent 

16016.32 

16000 

16.32 

Descent 

15497.29 

15500 

2.71 

Descent 

14992.09 

15000 

7.91 

Descent 

14499.40 

14500 

0.60 

Descent 

14023.03 

14000 

23.03 

Descent 

13521.65 

13500 

21.65 

Descent 

13019.29 

13000 

19.29 

Descent 

12515.42 

12500 

15.42 

Descent 

12020.22 

12000 

20.22 

Descent 

11521.54 

11500 

21.54 

Descent 

10975.08 

11000 

24.92 

Descent 

10516.08 

10500 

16.08 

Descent 

10012.22 

10000 

12.22 

Descent 

9503.27 

9500 

3.27 

Descent 

8989.56 

9000 

10.44 

Descent 

8503.66 

8500 

3.66 

Descent 

7996.16 

8000 

3.84 

Descent 

7505.77 

7500 

5.77 

Descent 

6990.75 

7000 

9.25 

Descent 

6493.86 

6500 

6.14 

Descent 

6006.63 

6000 

6.63 

Descent 

5494.39 

5500 

5.61 

Descent 

5009.86 

5000 

9.86 

Descent 

4502.20 

4500 

2.20 

Ascent 

2485.44 

2500 

14.56 

Aseent 

2989.83 

3000 

10.17 

Ascent 

3492.68 

3500 

7.32 

Ascent 

4016.26 

4000 

16.26 
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2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


Ascent 

4497.97 

4500 

2.03 

Ascent 

S028.65 

5000 

28.65 

Ascent 

5482.95 

5500 

17.05 

Ascent 

6009.95 

6000 

9.95 

Ascent 

6500.60 

6500 

0.60 

Ascent 

6977.07 

7000 

22.93 

Ascent 

7527.22 

7500 

27.22 

Ascent 

7978.51 

8000 

21.49 

Ascent 

8498.88 

8500 

1.12 

Ascent 

9009.58 

9000 

9.58 

Ascent 

9520.53 

9500 

20.53 

Ascent 

10026.01 

10000 

26.01 

Ascent 

10503.98 

10500 

3.98 

Ascent 

10998.33 

11000 

1.67 

Ascent 

11474.93 

11500 

25.07 

Ascent 

11983.53 

12000 

16.47 

Ascent 

12517.46 

12500 

17.46 

Ascent 

13015.16 

13000 

15.16 

Ascent 

13502.74 

13500 

2.74 

Ascent 

13996.69 

14000 

3.31 

Ascent 

14518.93 

14500 

18.93 

Ascent 

14982.53 

15000 

17.47 

Ascent 

15506.25 

15500 

6.25 

Ascent 

15979.88 

16000 

20.12 

Ascent 

16517.32 

16500 

17.32 

Ascent 

16985.93 

17000 

14.07 

Ascent 

17522.55 

17500 

22.55 

Ascent 

18021.79 

18000 

21.79 

Ascent 

18477.52 

18500 

22.48 

Ascent 

18982.21 

19000 

17.79 

Ascent 

19499.04 

19500 

0.96 

Ascent 

19989.47 

20000 

10.53 

Ascent 

20485.54 

20500 

14.46 

Ascent 

21006.30 

21000 

6.30 

Ascent 

21508.17 

21500 

8.17 

Ascent 

22011.40 

22000 

11.40 

Ascent 

22487.58 

22500 

12.42 

Ascent 

22984.24 

23000 

15.76 

Ascent 

23518.99 

23500 

18.99 

Ascent 

24001.08 

24000 

1.08 

Ascent 

24520.87 

24500 

20.87 

Ascent 

24982.20 

25000 

17.80 

Ascent 

25483.69 

25500 

16.31 

Ascent 

26000.69 

26000 

0.69 

Ascent 

26500.16 

26500 

0.16 
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K)  K)  K) 


2  Ascent  27005 

2  Ascent  27509 

2  Aseent  27981 

2  Ascent  28495 

2  Aseent  29002 

2  Ascent  29503 

Aseent  30006 

Aseent  30502 

Aseent  30993 

2  Ascent  31484 

2  Aseent  32005 

2  Aseent  32498 

2  Aseent  32991 

2  Aseent  33496 

2  Ascent  33998 

2  Aseent  34506 

2  Ascent  34999 

2  Aseent  35486 

2  Aseent  35993 

2  Ascent  36503 

2  Aseent  36989 

2  Aseent  37508 

2  Aseent  38007 

2  Aseent  38504 

2  Aseent  38996 

2  Aseent  39496 

2  Aseent  39997 

2  Aseent  40493 

2  Aseent  41004 

2  Aseent  41508 

2  Ascent  42005 

2  Aseent  42503 

2  Aseent  43001 

2  Aseent  43504 

2  Aseent  43992 

2  Ascent  44494 

2  Aseent  45004 

2  Aseent  45504 

2  Ascent  45992 

2  Aseent  46497 

2  Aseent  46996 

2  Ascent  47501 

2  Aseent  48000 

2  Ascent  48506 

2  Aseent  49004 


.39 

27000 

5.39 

.81 

27500 

9.81 

.07 

28000 

18.93 

.59 

28500 

4.41 

.04 

29000 

2.04 

.48 

29500 

3.48 

.91 

30000 

6.91 

.23 

30500 

2.23 

.97 

31000 

6.03 

.43 

31500 

15.57 

.32 

32000 

5.32 

.60 

32500 

1.40 

.54 

33000 

8.46 

.66 

33500 

3.34 

.75 

34000 

1.25 

.27 

34500 

6.27 

.97 

35000 

0.03 

.53 

35500 

13.47 

.27 

36000 

6.73 

.06 

36500 

3.06 

.82 

37000 

10.18 

.25 

37500 

8.25 

.46 

38000 

7.46 

.95 

38500 

4.95 

.73 

39000 

3.27 

.74 

39500 

3.26 

.79 

40000 

2.21 

.89 

40500 

6.11 

.08 

41000 

4.08 

.91 

41500 

8.91 

.59 

42000 

5.59 

.83 

42500 

3.83 

.27 

43000 

1.27 

.22 

43500 

4.22 

.31 

44000 

7.69 

.45 

44500 

5.55 

.24 

45000 

4.24 

.79 

45500 

4.79 

.61 

46000 

7.39 

.27 

46500 

2.73 

.07 

47000 

3.93 

.73 

47500 

1.73 

.67 

48000 

0.67 

.21 

48500 

6.21 

.06 

49000 

4.06 
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2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


Ascent 

49502.25 

49500 

2.25 

Ascent 

50006.58 

50000 

6.58 

Ascent 

50495.42 

50500 

4.58 

Ascent 

50996.03 

51000 

3.97 

Aseent 

51495.92 

51500 

4.08 

Ascent 

51994.73 

52000 

5.27 

Ascent 

52498.96 

52500 

1.04 

Aseent 

53001.66 

53000 

1.66 

Ascent 

53502.43 

53500 

2.43 

Ascent 

53997.13 

54000 

2.87 

Aseent 

54496.34 

54500 

3.66 

Ascent 

55000.10 

55000 

0.10 

Ascent 

55500.52 

55500 

0.52 

Ascent 

56001.12 

56000 

1.12 

Aseent 

56501.59 

56500 

1.59 

Aseent 

57001.64 

57000 

1.64 

Ascent 

57500.93 

57500 

0.93 

Ascent 

57999.11 

58000 

0.89 

Aseent 

58500.37 

58500 

0.37 

Aseent 

59004.63 

59000 

4.63 

Ascent 

59502.23 

59500 

2.23 

Ascent 

60002.20 

60000 

2.20 

Descent 

60002.20 

60000 

2.20 

Ascent 

2483.95 

2500 

16.05 

Ascent 

3024.69 

3000 

24.69 

Ascent 

3479.86 

3500 

20.14 

Aseent 

3977.73 

4000 

22.27 

Ascent 

4514.88 

4500 

14.88 

Ascent 

4985.18 

5000 

14.82 

Ascent 

5502.02 

5500 

2.02 

Ascent 

5991.15 

6000 

8.85 

Aseent 

6506.22 

6500 

6.22 

Ascent 

7003.30 

7000 

3.30 

Aseent 

7527.80 

7500 

27.80 

Aseent 

7986.75 

8000 

13.25 

Aseent 

8504.85 

8500 

4.85 

Ascent 

9005.94 

9000 

5.94 

Ascent 

9505.12 

9500 

5.12 

Ascent 

10000.32 

10000 

0.32 

Aseent 

10495.71 

10500 

4.29 

Ascent 

10986.70 

11000 

13.30 

Ascent 

11498.56 

11500 

1.44 

Aseent 

12004.87 

12000 

4.87 

Ascent 

12482.21 

12500 

17.79 

Ascent 

13000.69 

13000 

0.69 
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3 

Ascent 

13524.46 

13500 

24.46 

3 

Ascent 

13974.66 

14000 

25.34 

3 

Aseent 

14489.28 

14500 

10.72 

3 

Ascent 

15013.40 

15000 

13.40 

3 

Aseent 

15476.39 

15500 

23.61 

3 

Ascent 

15995.06 

16000 

4.94 

3 

Aseent 

16508.06 

16500 

8.06 

3 

Aseent 

17000.83 

17000 

0.83 

3 

Aseent 

17478.66 

17500 

21.34 

3 

Ascent 

17976.36 

18000 

23.64 

3 

Aseent 

18489.07 

18500 

10.93 

3 

Aseent 

18993.12 

19000 

6.88 

3 

Aseent 

19502.45 

19500 

2.45 

3 

Aseent 

19979.92 

20000 

20.08 

3 

Ascent 

20512.90 

20500 

12.90 

3 

Aseent 

21017.98 

21000 

17.98 

3 

Ascent 

21517.31 

21500 

17.31 

3 

Aseent 

21981.66 

22000 

18.34 

3 

Aseent 

22483.80 

22500 

16.20 

3 

Ascent 

23014.06 

23000 

14.06 

3 

Aseent 

23510.17 

23500 

10.17 

3 

Aseent 

24014.03 

24000 

14.03 

3 

Aseent 

24515.79 

24500 

15.79 

3 

Aseent 

25009.02 

25000 

9.02 

3 

Aseent 

25518.35 

25500 

18.35 

3 

Aseent 

26018.87 

26000 

18.87 

3 

Aseent 

26509.95 

26500 

9.95 

3 

Aseent 

26992.09 

27000 

7.91 

3 

Aseent 

27491.76 

27500 

8.24 

3 

Aseent 

27998.30 

28000 

1.70 

3 

Ascent 

28503.77 

28500 

3.77 

3 

Aseent 

29013.96 

29000 

13.96 

3 

Aseent 

29501.06 

29500 

1.06 

3 

Aseent 

29982.20 

30000 

17.80 

3 

Aseent 

30508.52 

30500 

8.52 

3 

Ascent 

31000.38 

31000 

0.38 

3 

Aseent 

31505.32 

31500 

5.32 

3 

Aseent 

31990.68 

32000 

9.32 

3 

Ascent 

32483.68 

32500 

16.32 

3 

Aseent 

33006.74 

33000 

6.74 

3 

Aseent 

33489.61 

33500 

10.39 

3 

Ascent 

33985.83 

34000 

14.17 

3 

Aseent 

34509.20 

34500 

9.20 

3 

Ascent 

35002.96 

35000 

2.96 

3 

Aseent 

35494.13 

35500 

5.87 
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3 

Ascent 

36007.22 

36000 

7.22 

3 

Ascent 

36493.53 

36500 

6.47 

3 

Aseent 

37004.46 

37000 

4.46 

3 

Ascent 

37514.91 

37500 

14.91 

3 

Aseent 

38000.63 

38000 

0.63 

3 

Ascent 

38490.97 

38500 

9.03 

3 

Aseent 

39002.10 

39000 

2.10 

3 

Aseent 

39505.91 

39500 

5.91 

3 

Aseent 

39999.67 

40000 

0.33 

3 

Ascent 

40509.29 

40500 

9.29 

3 

Aseent 

40994.23 

41000 

5.77 

3 

Aseent 

41494.77 

41500 

5.23 

3 

Aseent 

42003.53 

42000 

3.53 

3 

Aseent 

42501.71 

42500 

1.71 

3 

Ascent 

42996.93 

43000 

3.07 

3 

Aseent 

43490.89 

43500 

9.11 

3 

Ascent 

44008.25 

44000 

8.25 

3 

Aseent 

44499.12 

44500 

0.88 

3 

Aseent 

45004.24 

45000 

4.24 

3 

Ascent 

45499.89 

45500 

0.11 

3 

Aseent 

46000.13 

46000 

0.13 

3 

Aseent 

46502.41 

46500 

2.41 

3 

Aseent 

46996.07 

47000 

3.93 

3 

Aseent 

47504.43 

47500 

4.43 

3 

Aseent 

48000.67 

48000 

0.67 

3 

Aseent 

48500.55 

48500 

0.55 

3 

Aseent 

48998.27 

49000 

1.73 

3 

Aseent 

49505.21 

49500 

5.21 

3 

Aseent 

49997.46 

50000 

2.54 

3 

Aseent 

50501.64 

50500 

1.64 

3 

Ascent 

50999.21 

51000 

0.79 

3 

Aseent 

51499.18 

51500 

0.82 

3 

Aseent 

51998.08 

52000 

1.92 

3 

Aseent 

52502.38 

52500 

2.38 

3 

Aseent 

53008.68 

53000 

8.68 

3 

Ascent 

53498.83 

53500 

1.17 

3 

Aseent 

53997.13 

54000 

2.87 

3 

Aseent 

54503.88 

54500 

3.88 

3 

Ascent 

54996.24 

55000 

3.76 

3 

Aseent 

55500.52 

55500 

0.52 

3 

Aseent 

56001.12 

56000 

1.12 

3 

Ascent 

56501.59 

56500 

1.59 

3 

Aseent 

57001.64 

57000 

1.64 

3 

Ascent 

57500.93 

57500 

0.93 

3 

Aseent 

57999.11 

58000 

0.89 
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3 

Ascent 

58500.37 

58500 

0.37 

3 

Ascent 

58999.95 

59000 

0.05 

3 

Aseent 

59502.23 

59500 

2.23 

3 

Ascent 

60002.20 

60000 

2.20 

3 

Descent 

60002.20 

60000 

2.20 

3 

Descent 

59502.23 

59500 

2.23 

3 

Descent 

59004.63 

59000 

4.63 

3 

Descent 

58495.80 

58500 

4.20 

3 

Descent 

57999.11 

58000 

0.89 

3 

Descent 

57487.87 

57500 

12.13 

3 

Descent 

57001.64 

57000 

1.64 

3 

Descent 

56509.90 

56500 

9.90 

3 

Descent 

56013.28 

56000 

13.28 

3 

Descent 

55492.61 

55500 

7.39 

3 

Descent 

54996.24 

55000 

3.76 

3 

Descent 

54477.50 

54500 

22.50 

3 

Descent 

53997.13 

54000 

2.87 

3 

Descent 

53473.69 

53500 

26.31 

3 

Deseent 

52987.63 

53000 

12.37 

3 

Descent 

52526.38 

52500 

26.38 

3 

Descent 

52004.77 

52000 

4.77 

3 

Descent 

51518.78 

51500 

18.78 

3 

Descent 

51015.16 

51000 

15.16 

3 

Deseent 

50520.32 

50500 

20.32 

3 

Descent 

50018.73 

50000 

18.73 

3 

Descent 

49514.11 

49500 

14.11 

3 

Deseent 

49009.85 

49000 

9.85 

3 

Descent 

48503.38 

48500 

3.38 

3 

Descent 

47973.09 

48000 

26.91 

3 

Descent 

47507.13 

47500 

7.13 

3 

Descent 

46969.80 

47000 

30.20 

3 

Descent 

46512.68 

46500 

12.68 

3 

Descent 

45975.07 

46000 

24.93 

3 

Descent 

45473.00 

45500 

27.00 

3 

Descent 

45016.19 

45000 

16.19 

3 

Descent 

44499.12 

44500 

0.88 

3 

Descent 

44017.36 

44000 

17.36 

3 

Descent 

43510.89 

43500 

10.89 

3 

Descent 

43005.61 

43000 

5.61 

3 

Descent 

42520.79 

42500 

20.79 

3 

Descent 

42018.01 

42000 

18.01 

3 

Descent 

41502.85 

41500 

2.85 

3 

Descent 

41004.08 

41000 

4.08 

3 

Descent 

40490.04 

40500 

9.96 

3 

Descent 

39982.77 

40000 

17.23 
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3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


Descent 

39476.58 

39500 

23.42 

Descent 

39005.68 

39000 

5.68 

Descent 

38504.95 

38500 

4.95 

Descent 

38000.63 

38000 

0.63 

Descent 

37499.91 

37500 

0.09 

Descent 

37015.85 

37000 

15.85 

Descent 

36498.29 

36500 

1.71 

Descent 

35999.46 

36000 

0.54 

Descent 

35489.56 

35500 

10.44 

Descent 

34999.97 

35000 

0.03 

Descent 

34507.73 

34500 

7.73 

Descent 

34024.62 

34000 

24.62 

Descent 

33503.71 

33500 

3.71 

Descent 

32977.72 

33000 

22.28 

Descent 

32493.17 

32500 

6.83 

Descent 

31989.35 

32000 

10.65 

Descent 

31477.89 

31500 

22.11 

Descent 

31014.50 

31000 

14.50 

Descent 

30502.23 

30500 

2.23 

Descent 

29988.38 

30000 

11.62 

Descent 

29496.20 

29500 

3.80 

Descent 

28978.23 

29000 

21.77 

Descent 

28500.26 

28500 

0.26 

Descent 

27969.59 

28000 

30.41 

Descent 

27522.23 

27500 

22.23 

Descent 

26978.80 

27000 

21.20 

Descent 

26478.39 

26500 

21.61 

Descent 

25971.83 

26000 

28.17 

Descent 

25471.09 

25500 

28.91 

Descent 

24977.04 

25000 

22.96 

Descent 

24489.43 

24500 

10.57 

Descent 

23992.12 

24000 

7.88 

Descent 

23509.19 

23500 

9.19 

Descent 

23026.58 

23000 

26.58 

Descent 

22504.60 

22500 

4.60 

Descent 

21982.59 

22000 

17.41 

Descent 

21499.94 

21500 

0.06 

Descent 

21021.58 

21000 

21.58 

Descent 

20492.60 

20500 

7.40 

Descent 

19987.73 

20000 

12.27 

Descent 

19505.02 

19500 

5.02 

Descent 

19001.51 

19000 

1.51 

Descent 

18489.89 

18500 

10.11 

Descent 

17922.89 

17900 

22.89 

Descent 

17521.75 

17500 

21.75 
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3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


Descent 

17006.33 

17000 

6.33 

Descent 

16509.60 

16500 

9.60 

Descent 

15973.06 

16000 

26.94 

Descent 

15494.31 

15500 

5.69 

Descent 

15014.13 

15000 

14.13 

Descent 

14508.08 

14500 

8.08 

Descent 

14000.96 

14000 

0.96 

Descent 

13478.95 

13500 

21.05 

Descent 

12978.65 

13000 

21.35 

Descent 

12480.86 

12500 

19.14 

Descent 

11988.20 

12000 

11.80 

Descent 

11502.50 

11500 

2.50 

Descent 

10982.83 

11000 

17.17 

Descent 

10474.73 

10500 

25.27 

Descent 

10015.35 

10000 

15.35 

Descent 

9484.16 

9500 

15.84 

Descent 

9004.73 

9000 

4.73 

Descent 

8509.64 

8500 

9.64 

Descent 

8012.64 

8000 

12.64 

Descent 

7504.04 

7500 

4.04 

Descent 

6998.17 

7000 

1.83 

Descent 

6492.17 

6500 

7.83 

Descent 

5995.01 

6000 

4.99 

Descent 

5486.22 

5500 

13.78 

Descent 

4991.07 

5000 

8.93 

Descent 

4508.01 

4500 

8.01 

Ascent 

3013.57 

3000 

13.57 

Ascent 

3502.94 

3500 

2.94 

Ascent 

3979.81 

4000 

20.19 

Ascent 

4512.24 

4500 

12.24 

Ascent 

5012.01 

5000 

12.01 

Ascent 

5487.31 

5500 

12.69 

Ascent 

6011.06 

6000 

11.06 

Ascent 

6496.67 

6500 

3.33 

Ascent 

6978.20 

7000 

21.80 

Ascent 

7505.77 

7500 

5.77 

Ascent 

7989.69 

8000 

10.31 

Ascent 

8476.78 

8500 

23.22 

Ascent 

8994.41 

9000 

5.59 

Ascent 

9502.03 

9500 

2.03 

Ascent 

9980.29 

10000 

19.71 

Ascent 

10483.63 

10500 

16.37 

Ascent 

10986.05 

11000 

13.95 

Ascent 

11495.94 

11500 

4.06 

Ascent 

12028.23 

12000 

28.23 
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4 

Ascent 

12526.28 

12500 

26.28 

4 

Ascent 

12998.63 

13000 

1.37 

4 

Aseent 

13483.85 

13500 

16.15 

4 

Ascent 

14015.91 

14000 

15.91 

4 

Aseent 

14510.25 

14500 

10.25 

4 

Ascent 

14984.00 

15000 

16.00 

4 

Aseent 

15512.24 

15500 

12.24 

4 

Aseent 

16007.21 

16000 

7.21 

4 

Aseent 

16502.65 

16500 

2.65 

4 

Ascent 

17003.19 

17000 

3.19 

4 

Aseent 

17516.96 

17500 

16.96 

4 

Aseent 

18017.73 

18000 

17.73 

4 

Aseent 

18516.30 

18500 

16.30 

4 

Aseent 

18996.47 

19000 

3.53 

4 

Ascent 

19493.07 

19500 

6.93 

4 

Aseent 

19992.94 

20000 

7.06 

4 

Ascent 

20507.61 

20500 

7.61 

4 

Aseent 

20990.14 

21000 

9.86 

4 

Aseent 

21491.72 

21500 

8.28 

4 

Ascent 

22010.47 

22000 

10.47 

4 

Aseent 

22492.30 

22500 

7.70 

4 

Aseent 

22982.31 

23000 

17.69 

4 

Aseent 

23494.50 

23500 

5.50 

4 

Aseent 

23992.12 

24000 

7.88 

4 

Aseent 

24509.71 

24500 

9.71 

4 

Aseent 

25010.05 

25000 

10.05 

4 

Aseent 

25499.43 

25500 

0.57 

4 

Aseent 

26017.80 

26000 

17.80 

4 

Aseent 

26493.62 

26500 

6.38 

4 

Aseent 

26989.88 

27000 

10.12 

4 

Ascent 

27506.43 

27500 

6.43 

4 

Aseent 

28000.60 

28000 

0.60 

4 

Aseent 

28503.77 

28500 

3.77 

4 

Aseent 

28998.47 

29000 

1.53 

4 

Aseent 

29494.99 

29500 

5.01 

4 

Ascent 

29995.79 

30000 

4.21 

4 

Aseent 

30504.75 

30500 

4.75 

4 

Aseent 

31001.67 

31000 

1.67 

4 

Ascent 

31489.64 

31500 

10.36 

4 

Aseent 

32010.64 

32000 

10.64 

4 

Aseent 

32510.80 

32500 

10.80 

4 

Ascent 

32992.91 

33000 

7.09 

4 

Aseent 

33492.43 

33500 

7.57 

4 

Ascent 

34008.80 

34000 

8.80 

4 

Aseent 

34498.95 

34500 

1.05 

91 


4 

Ascent 

34991.02 

35000 

8.98 

4 

Ascent 

35507.82 

35500 

7.82 

4 

Aseent 

36001.02 

36000 

1.02 

4 

Ascent 

36495.12 

36500 

4.88 

4 

Aseent 

36991.45 

37000 

8.55 

4 

Ascent 

37498.25 

37500 

1.75 

4 

Aseent 

38000.63 

38000 

0.63 

4 

Aseent 

38503.20 

38500 

3.20 

4 

Aseent 

39002.10 

39000 

2.10 

4 

Ascent 

39491.24 

39500 

8.76 

4 

Aseent 

39999.67 

40000 

0.33 

4 

Aseent 

40497.74 

40500 

2.26 

4 

Aseent 

40996.20 

41000 

3.80 

4 

Aseent 

41502.85 

41500 

2.85 

4 

Ascent 

42001.46 

42000 

1.46 

4 

Aseent 

42505.95 

42500 

5.95 

4 

Ascent 

43001.27 

43000 

1.27 

4 

Aseent 

43504.22 

43500 

4.22 

4 

Aseent 

44001.41 

44000 

1.41 

4 

Ascent 

44503.78 

44500 

3.78 

4 

Aseent 

45001.85 

45000 

1.85 

4 

Aseent 

45499.89 

45500 

0.11 

4 

Aseent 

45997.62 

46000 

2.38 

4 

Aseent 

46497.27 

46500 

2.73 

4 

Aseent 

47001.33 

47000 

1.33 

4 

Aseent 

47499.04 

47500 

0.96 

4 

Aseent 

47997.91 

48000 

2.09 

4 

Aseent 

48497.73 

48500 

2.27 

4 

Aseent 

49004.06 

49000 

4.06 

4 

Aseent 

49496.32 

49500 

3.68 

4 

Ascent 

49997.46 

50000 

2.54 

4 

Aseent 

50501.64 

50500 

1.64 

4 

Aseent 

50999.21 

51000 

0.79 

4 

Aseent 

51502.45 

51500 

2.45 

4 

Aseent 

52001.43 

52000 

1.43 

4 

Ascent 

52495.53 

52500 

4.47 

4 

Aseent 

52998.15 

53000 

1.85 

4 

Aseent 

53498.83 

53500 

1.17 

4 

Ascent 

53997.13 

54000 

2.87 

4 

Aseent 

54500.11 

54500 

0.11 

4 

Aseent 

55000.10 

55000 

0.10 

4 

Ascent 

55500.52 

55500 

0.52 

4 

Aseent 

55997.06 

56000 

2.94 

4 

Aseent 

56497.44 

56500 

2.56 

4 

Descent 

56501.59 

56500 

1.59 
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4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 


Descent 

55997.06 

56000 

2.94 

Descent 

55508.44 

55500 

8.44 

Descent 

54996.24 

55000 

3.76 

Descent 

54500.11 

54500 

0.11 

Descent 

54000.81 

54000 

0.81 

Descent 

53491.64 

53500 

8.36 

Descent 

53008.68 

53000 

8.68 

Descent 

52485.26 

52500 

14.74 

Descent 

51994.73 

52000 

5.27 

Descent 

51495.92 

51500 

4.08 

Descent 

51008.78 

51000 

8.78 

Descent 

S0479.88 

50500 

20.13 

Descent 

49982.28 

50000 

17.72 

Descent 

49511.14 

49500 

11.14 

Descent 

48989.59 

49000 

10.41 

Descent 

48500.55 

48500 

0.55 

Descent 

48017.23 

48000 

17.23 

Ascent 

2514.32 

2500 

14.32 

Ascent 

2986.79 

3000 

13.21 

Ascent 

3496.27 

3500 

3.73 

Ascent 

4519.11 

4500 

19.11 

Ascent 

5005.57 

5000 

5.57 

Ascent 

5503.65 

5500 

3.65 

Ascent 

5989.48 

6000 

10.52 

Ascent 

6509.03 

6500 

9.03 

Ascent 

7013.57 

7000 

13.57 

Ascent 

7488.39 

7500 

11.61 

Ascent 

7986.75 

8000 

13.25 

Ascent 

8493.51 

8500 

6.49 

Ascent 

8979.25 

9000 

20.75 

Ascent 

9511.90 

9500 

11.90 

Ascent 

9981.54 

10000 

18.46 

Ascent 

10498.89 

10500 

1.11 

Ascent 

10995.74 

11000 

4.26 

Ascent 

11518.92 

11500 

18.92 

Ascent 

11985.53 

12000 

14.47 

Ascent 

12525.60 

12500 

25.60 

Ascent 

12996.55 

13000 

3.45 

Ascent 

13497.15 

13500 

2.85 

Ascent 

14028.72 

14000 

28.72 

Ascent 

14481.33 

14500 

18.67 

Ascent 

14994.29 

15000 

5.71 

Ascent 

15504.02 

15500 

4.02 

Ascent 

15990.51 

16000 

9.49 

Ascent 

16505.74 

16500 

5.74 
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5 

Ascent 

17000.83 

17000 

0.83 

5 

Ascent 

17482.65 

17500 

17.35 

5 

Aseent 

17988.52 

18000 

11.48 

5 

Ascent 

18485.77 

18500 

14.23 

5 

Aseent 

19008.22 

19000 

8.22 

5 

Ascent 

19500.74 

19500 

0.74 

5 

Aseent 

19985.13 

20000 

14.88 

5 

Aseent 

20521.74 

20500 

21.74 

5 

Aseent 

21008.10 

21000 

8.10 

5 

Ascent 

21510.91 

21500 

10.91 

5 

Aseent 

22004.90 

22000 

4.90 

5 

Aseent 

22499.86 

22500 

0.14 

5 

Aseent 

23011.17 

23000 

11.17 

5 

Aseent 

23508.21 

23500 

8.21 

5 

Ascent 

24013.04 

24000 

13.04 

5 

Aseent 

24479.30 

24500 

20.70 

5 

Ascent 

25011.08 

25000 

11.08 

5 

Aseent 

25485.79 

25500 

14.21 

5 

Aseent 

26004.96 

26000 

4.96 

5 

Ascent 

26499.07 

26500 

0.93 

5 

Aseent 

26995.41 

27000 

4.59 

5 

Aseent 

27497.39 

27500 

2.61 

5 

Aseent 

28012.08 

28000 

12.08 

5 

Aseent 

28483.89 

28500 

16.11 

5 

Aseent 

29010.39 

29000 

10.39 

5 

Aseent 

29494.99 

29500 

5.01 

5 

Aseent 

29993.32 

30000 

6.68 

5 

Aseent 

30492.15 

30500 

7.85 

5 

Aseent 

30991.40 

31000 

8.60 

5 

Aseent 

31502.71 

31500 

2.71 

5 

Ascent 

31998.66 

32000 

1.34 

5 

Aseent 

32499.95 

32500 

0.05 

5 

Aseent 

33012.27 

33000 

12.27 

5 

Aseent 

33500.89 

33500 

0.89 

5 

Aseent 

33998.75 

34000 

1.25 

5 

Ascent 

34509.20 

34500 

9.20 

5 

Aseent 

35010.42 

35000 

10.42 

5 

Aseent 

35497.17 

35500 

2.83 

5 

Ascent 

35994.82 

36000 

5.18 

5 

Aseent 

36501.46 

36500 

1.46 

5 

Aseent 

37007.71 

37000 

7.71 

5 

Ascent 

37498.25 

37500 

1.75 

5 

Aseent 

37988.69 

38000 

11.31 

5 

Ascent 

38490.97 

38500 

9.03 

5 

Aseent 

38993.14 

39000 

6.86 
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5 

Ascent 

39502.24 

39500 

2.24 

5 

Ascent 

39997.79 

40000 

2.21 

5 

Aseent 

40509.29 

40500 

9.29 

5 

Ascent 

41004.08 

41000 

4.08 

5 

Aseent 

41506.89 

41500 

6.89 

5 

Ascent 

42007.66 

42000 

7.66 

5 

Aseent 

42505.95 

42500 

5.95 

5 

Aseent 

42992.59 

43000 

7.41 

5 

Aseent 

43504.22 

43500 

4.22 

5 

Ascent 

43999.14 

44000 

0.86 

5 

Aseent 

44503.78 

44500 

3.78 

5 

Aseent 

44999.46 

45000 

0.54 

5 

Aseent 

45499.89 

45500 

0.11 

5 

Aseent 

45997.62 

46000 

2.38 

5 

Ascent 

46499.84 

46500 

0.16 

5 

Aseent 

47003.96 

47000 

3.96 

5 

Ascent 

47496.35 

47500 

3.65 

5 

Aseent 

47997.91 

48000 

2.09 

5 

Aseent 

48506.21 

48500 

6.21 

5 

Ascent 

48998.27 

49000 

1.73 

5 

Aseent 

49496.32 

49500 

3.68 

5 

Aseent 

49997.46 

50000 

2.54 

5 

Aseent 

50498.54 

50500 

1.46 

5 

Aseent 

50999.21 

51000 

0.79 

5 

Aseent 

51499.18 

51500 

0.82 

5 

Aseent 

51998.08 

52000 

1.92 

5 

Aseent 

52502.38 

52500 

2.38 

5 

Aseent 

53001.66 

53000 

1.66 

5 

Aseent 

53498.83 

53500 

1.17 

5 

Aseent 

53997.13 

54000 

2.87 

5 

Ascent 

54496.34 

54500 

3.66 

5 

Aseent 

55000.10 

55000 

0.10 

5 

Aseent 

55500.52 

55500 

0.52 

5 

Aseent 

56001.12 

56000 

1.12 

5 

Aseent 

56501.59 

56500 

1.59 

5 

Ascent 

57001.64 

57000 

1.64 

5 

Aseent 

57500.93 

57500 

0.93 

5 

Aseent 

57994.65 

58000 

5.35 

5 

Ascent 

58500.37 

58500 

0.37 

5 

Aseent 

58999.95 

59000 

0.05 

5 

Aseent 

59502.23 

59500 

2.23 

5 

Ascent 

60002.20 

60000 

2.20 

5 

Aseent 

60499.37 

60500 

0.63 

5 

Ascent 

61003.55 

61000 

3.55 

5 

Aseent 

61509.70 

61500 

9.70 
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5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 


Ascent 

62001.41 

62000 

1.41 

Descent 

62001.41 

62000 

1.41 

Descent 

61499.14 

61500 

0.86 

Descent 

60998.40 

61000 

1.60 

Descent 

60504.40 

60500 

4.40 

Descent 

60002.20 

60000 

2.20 

Descent 

59497.44 

59500 

2.56 

Descent 

59014.00 

59000 

14.00 

Descent 

58509.51 

58500 

9.51 

Descent 

57985.73 

58000 

14.27 

Descent 

57483.52 

57500 

16.48 

Descent 

56980.39 

57000 

19.61 

Descent 

56489.14 

56500 

10.86 

Descent 

56017.34 

56000 

17.34 

Descent 

55468.89 

55500 

31.11 

Descent 

S5038.77 

55000 

38.771 

Descent 

54515.20 

54500 

15.20 

Descent 

54030.29 

54000 

30.29 

Descent 

53520.41 

53500 

20.41 

Descent 

53033.27 

53000 

33.27 

Descent 

52516.09 

52500 

16.09 

Descent 

52004.77 

52000 

4.77 

Descent 

51486.13 

51500 

13.87 

Descent 

51034.30 

51000 

34.30 

Descent 

50529.67 

50500 

29.67 

Descent 

49973.18 

50000 

26.82 

Descent 

49472.61 

49500 

27.39 

Descent 

48969.34 

49000 

30.66 

Descent 

48520.35 

48500 

20.35 

Descent 

48006.19 

48000 

6.19 

Descent 

47515.21 

47500 

15.21 

Descent 

47001.33 

47000 

1.33 

Descent 

46479.31 

46500 

20.69 

Descent 

46005.14 

46000 

5.14 

Descent 

45536.63 

45500 

36.63 

Descent 

45030.54 

45000 

30.54 

Descent 

44468.82 

44500 

31.18 

Descent 

43980.93 

44000 

19.07 

Descent 

43519.79 

43500 

19.79 

Descent 

42966.58 

43000 

33.42 

Descent 

42497.48 

42500 

2.52 

Descent 

41974.59 

42000 

25.41 

Descent 

41478.63 

41500 

21.37 

Descent 

41025.77 

41000 

25.77 

Descent 

40501.59 

40500 

1.59 
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5 

Descent 

40022.22 

40000 

22.22 

5 

Descent 

39520.59 

39500 

20.59 

5 

Descent 

38984.20 

39000 

15.80 

5 

Descent 

38483.98 

38500 

16.02 

5 

Descent 

37968.23 

38000 

31.77 

5 

Descent 

37508.25 

37500 

8.25 

5 

Descent 

37019.11 

37000 

19.11 

5 

Descent 

36504.64 

36500 

4.64 

5 

Descent 

35996.37 

36000 

3.63 

5 

Descent 

35474.36 

35500 

25.64 

5 

Descent 

34996.98 

35000 

3.02 

5 

Descent 

34487.24 

34500 

12.76 

5 

Descent 

33971.48 

34000 

28.52 

5 

Descent 

33516.39 

33500 

16.39 

5 

Descent 

33012.27 

33000 

12.27 

5 

Descent 

32498.60 

32500 

1.40 

5 

Descent 

32018.65 

32000 

18.65 

5 

Descent 

31506.63 

31500 

6.63 

5 

Descent 

30999.10 

31000 

0.90 

5 

Descent 

30511.04 

30500 

11.04 

5 

Descent 

29980.96 

30000 

19.04 

5 

Descent 

29511.99 

29500 

11.99 

5 

Descent 

28980.60 

29000 

19.40 

5 

Descent 

28500.27 

28500 

0.27 

5 

Descent 

27974.18 

28000 

25.82 

5 

Descent 

27501.91 

27500 

1.91 

5 

Descent 

27008.72 

27000 

8.72 

5 

Descent 

26505.60 

26500 

5.60 

5 

Descent 

26012.45 

26000 

12.45 

5 

Descent 

25484.74 

25500 

15.26 

5 

Descent 

25000.77 

25000 

0.77 

5 

Descent 

24509.71 

24500 

9.71 

5 

Descent 

24036.97 

24000 

36.97 

5 

Descent 

23504.30 

23500 

4.30 

5 

Descent 

23027.54 

23000 

27.54 

5 

Descent 

22477.19 

22500 

22.81 

5 

Descent 

22006.75 

22000 

6.75 

5 

Descent 

21520.05 

21500 

20.05 

5 

Descent 

21007.20 

21000 

7.20 

5 

Descent 

20491.72 

20500 

8.28 

5 

Descent 

20005.96 

20000 

5.96 

5 

Descent 

19495.63 

19500 

4.37 

5 

Descent 

19010.74 

19000 

10.74 

5 

Descent 

18497.31 

18500 

2.69 

5 

Descent 

17982.03 

18000 

17.97 
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5 

Descent 

17518.55 

17500 

18.55 

5 

Descent 

17019.68 

17000 

19.68 

5 

Descent 

16483.36 

16500 

16.64 

5 

Descent 

15988.99 

16000 

11.01 

5 

Descent 

15519.71 

15500 

19.71 

5 

Descent 

14988.41 

15000 

11.59 

5 

Descent 

14494.35 

14500 

5.65 

5 

Descent 

13988.87 

14000 

11.13 

5 

Descent 

13464.97 

13500 

35.03 

5 

Descent 

12984.17 

13000 

15.83 

5 

Descent 

12529.67 

12500 

29.67 

5 

Descent 

12006.21 

12000 

6.21 

5 

Descent 

11507.75 

11500 

7.75 

5 

Descent 

10980.24 

11000 

19.76 

5 

Descent 

10481.73 

10500 

18.27 

5 

Descent 

10029.76 

10000 

29.76 

5 

Descent 

9502.03 

9500 

2.03 

5 

Descent 

9012.01 

9000 

12.01 

5 

Descent 

8488.12 

8500 

11.88 

5 

Descent 

8009.10 

8000 

9.10 

5 

Descent 

7500.56 

7500 

0.56 

5 

Descent 

6999.30 

7000 

0.70 

5 

Descent 

6493.30 

6500 

6.70 

5 

Descent 

5992.25 

6000 

7.75 

5 

Descent 

5494.93 

5500 

5.07 

5 

Descent 

5009.86 

5000 

9.86 
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Appendix  D;  Global  Hawk  Data  Reports  &  Flight  Mode 


The  following  are  notes  pertaining  to  the  Global  Hawk  data  reports  based  on 
indieated  flight  mode  (Climb,  Altitude  Hold,  and  Deseend).  Approximate  Flight  Level  is 
reported  in  thousands  of  feet.  Aetual  Flight  Mode  indieates  any  possible  diserepaneies  in 
Global  Hawk  data  reports  based  on  differenees  from  Indieated  Flight  Mode.  Any  data 
reports  used  in  the  statistieal  analysis  of  this  researeh  are  indieated  in  the  last  eolumn. 


Mission 

# 

Indicated 

Flight 

Mode 

Time  (sec) 

Approximate 
Flight  Levels 
(Kft) 

Actual  Flight  Mode 

Data 

Used 

1 

Climb 

57898-58238 

2.3-19.7 

Yes 

1 

Alt  Hold 

58239-58617 

18.3-20.2 

Climb  &  Descend 

Yes 

1 

Deseend 

58631-58698 

18.3-15.5 

Yes 

1 

Alt  Hold 

58699-58769 

15.5-14.3 

Descend 

Yes 

1 

Deseend 

58770-58854 

14.3-10.2 

Yes 

1 

Alt  Hold 

58855-59500 

10.2-4.3 

Descend 

Yes 

2 

Climb 

60668-61254 

2.4-30.5 

Yes 

2 

Alt  Hold 

61255-61533 

30.5-31.0 

Climb  &  Descend 

Yes 

2 

Climb 

61534-61583 

31.0-32.5 

Yes 

2 

Alt  Hold 

61584-61702 

32.5-35.0 

Climb 

Yes 

2 

Alt  Hold 

61703-62016 

35.0-35.07 

Climb  &  Descend 

No 

2 

Climb 

62017-63559 

35.07-55.0 

Yes 

2 

Climb 

63560-64926 

55.0-59.0 

Slow;  Climb 

Yes 

2 

Climb 

64927-66791 

59.0-60.0 

Slow;  Climb  &  Descend 

Yes 

2 

Climb 

66792-72190 

59.7-60.4 

Climb  &  Descend 

No 

2 

Deseend 

72192-72344 

59.6-60.4 

Yes 

3 

Climb 

56236-56831 

2.1-30.5 

Yes 

3 

Alt  Hold 

56832-57137 

30.5-31.1 

Climb  &  Descend 

Yes 

3 

Climb 

57138-57250 

31.0-34.5 

Yes 

3 

Alt  Hold 

57251-57520 

34.5-35.0 

Climb  &  Descend 

Yes 

3 

Climb 

57521-59898 

35.0-58.0 

Yes 

3 

Climb 

59899-60260 

58.0-59.0 

Slow;  Climb 

Yes 

3 

Climb 

60261-66256 

59.0-60.0 

Slow;  Climb  &  Descend 

Yes 

3 

Climb 

66257-67446 

59.9-60.1 

Climb  &  Descend 

Yes 

3 

Deseend 

67447-68213 

37.7-59.9 

Yes 
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3 

Alt  Hold 

68214-68412 

37.5-37.7 

Climb  &  Descend 

Yes 

3 

Descend 

68413-68461 

35.2-37.5 

Yes 

3 

Alt  Hold 

68462-68577 

35.0-35.2 

Climb  &  Descend 

Yes 

3 

Descend 

68578-68700 

28.7-35.0 

Yes 

3 

Alt  Hold 

68701-68887 

28.5-28.7 

Climb  &  Descend 

Yes 

3 

Descend 

68888-68954 

24.7-28.5 

Yes 

3 

Alt  Hold 

68955-69054 

24.4-24.7 

Climb  &  Descend 

Yes 

3 

Descend 

69055-69127 

20.4-24.4 

Yes 

3 

Alt  Hold 

69128-69381 

18.3-20.3 

Climb  &  Descend 

Yes 

3 

Descend 

69382-69435 

15.3-18.3 

Yes 

3 

Alt  Hold 

69436-69535 

14.2-15.3 

Climb  &  Descend 

Yes 

3 

Descend 

69536-69614 

10.1-14.1 

Yes 

3 

Alt  Hold 

69615-69928 

4.3-10.0 

Climb  &  Descend 

Yes 

4 

Climb 

52017-55314 

2.2-54.0 

Yes 

4 

Climb 

55315-56585 

54.0-55.0 

Slow;  Climb  &  Descend 

Yes 

4 

Climb 

56586-56767 

54.9-55.0 

Climb  &  Descend 

No 

4 

Climb 

56768-58244 

55.0-56.0 

Slow;  Climb  &  Descend 

Yes 

4 

Climb 

58245-66218 

55.6-56.9 

Climb  &  Descend 

Yes 

4 

Descend 

66219-66258 

56.8-56.8 

Climb  &  Descend 

Yes 

4 

Descend 

66259-66614 

47.9-56.8 

Yes 

5 

Climb 

53210-56095 

2.2-57.0 

Yes 

5 

Climb 

56096-56489 

57.0-58.0 

Slow;  Climb  &  Descend 

Yes 

5 

Climb 

56490-58041 

58.0-59.0 

Slow;  Climb  &  Descend 

Yes 

5 

Climb 

58042-63062 

59.0-60.0 

Slow;  Climb  &  Descend 

Yes 

5 

Climb 

63063-85216 

59.7-62.3 

Slow;  Climb  &  Descend 

Yes 

5 

Descend 

85217-85634 

60.0-62.1 

Slow;  Climb  &  Descend 

Yes 

5 

Descend 

85635-86094 

35.2-60.0 

Climb  &  Descend 

Yes 

5 

Alt  Hold 

86095-86223 

35.0-35.2 

Climb  &  Descend 

No 

5 

Descend 

86224-86335 

28.7-35.0 

Yes 

5 

Alt  Hold 

86336-86488 

28.5-28.7 

Climb  &  Descend 

No 

5 

Descend 

86489-86551 

24.6-28.5 

Yes 

5 

Alt  Hold 

86552-86649 

24.4-24.6 

Climb  &  Descend 

Yes 

5 

Descend 

86650-86714 

20.4-24.4 

Yes 

5 

Alt  Hold 

86715-86834 

20.1-20.4 

Climb  &  Descend 

No 

5 

Alt  Hold 

86835-86904 

18.3-20.1 

Descend 

Yes 

5 

Alt  Hold 

86905-86999 

18.3-18.4 

Climb  &  Descend 

No 

5 

Descend 

87000-87049 

15.4-18.3 

Yes 

5 

Alt  Hold 

87050-87120 

15.2-15.4 

Climb  &  Descend 

No 

5 

Alt  Hold 

87121-87156 

14.2-15.2 

Descend 

Yes 

5 

Descend 

87157-87231 

10.0-14.2 

Yes 

5 

Alt  Hold 

87232-87316 

9.8-10.0 

Climb  &  Descend 

No 

5 

Alt  Hold 

87317-87573 

4.3-9.8 

Descend 

Yes 

100 


Appendix  E:  Global  Hawk  Dimensions 


The  following  diagram  illustrates  the  overall  dimensions  of  the  Global  Hawk 
airframe.  The  loeation  of  the  air  data  sensors  on  the  tails  of  the  aireraft  are  denoted 
(Northrop  Grumman  2001). 


11 6.2  ft 
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